AD-A091  073  CENTER  FOR  NAVAL  ANALYSES  ALEXANDRIA  VA  OPERATIONS  EV— ETC  F/6  1/3 

A  REPRINT  OF  *A  METHOD  OF  ESTIMATING  PLANE  VULNERABILITY  BASED  — ETC(U) 
JUL  80  A  WALD 
UNCLASSIFIED  CRC-432 


I 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURLAU  OF  STANDARDS  1963  A 


Abraham  Wald 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh*fi  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  0OVT  ACCESSION  NO. 

CRC  432  p' 

S.  RECIPIENT'S  CATALOG  NUMBER 

'  ^73' 

4.  TITLE  ranrf  SuHItl,) 

A  REPRINT  OF  "A  METHOD  OF  ESTIMATING  PLANE 
VULNERABILITY  BASED  ON  DAMAGE  OF  SURVIVORS" 

BY  ABRAHAM  WALD 

s.  TYPE  OF  REPORT  *  PERIOD  COVERED 

Reprint  -  1943 

1  PERFORMING  ORG.  REPORT  NUMBER 

See  Block  18 

7.  AUTHORO) 

j  Abraham  aid 

1.  CONTRACT  OR  GRANT  HUMBERTS 

S  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Statistical  Research  Group/National 

Defense  Research  Committee, 

Applied  Mathematics  Panel 

to.  program  element,  project,  task 

AREA  4  WORK  UNIT  NUMBERS 

1  1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Center  for  Naval  Analyses 

2000  N.  Beauregard  Street 

Alexandria,  Virginia  22311 

12.  REPORT  DATE 

juiy  1980 

13  BE*  OF  PAGES 

Q1 

•  4.  MONITORING  AGENCY  NAME  4  ADDRES$</<  dlllottnl  Iran  Controlling  Otllco) 

Office  of  the  Chief  of  Naval  Operations  (Op03) 
Department  of  the  Navy 

Washington,  D.C.  20350 

IS.  SECURITY  CLASS,  (ol  thlo  roport) 

Unclassified 

IS*.  DECL  ASSlFIC  ATION/ DOWNGRADING 
SCHEDULE 

1  1«.  OISTRIBUTiqN  STATEMENT  (of  thh  Report)  1 

Approved  for  public  release;  distribution  unlimited 


'7.  DISTRIBUTION  STATEMENT  fol  the  mbetrect  entered  In  Block  20,  II  dllterent  from  Report) 


IS.  SUPPLEMENTARY  NOTES 

1)  SRG  Memo  85/ AMP  Memo  76.1 

2)  SRG  Memo  87/ AMP  Memo  76.2 

3)  SRG  Memo  88/AMP  Memo  76.3 

4)  SRG  Memo  89/AMP  Memo  76.4 


5)  SRG  Memo  96/AMP  Memo  76.5 

6)  SRG  Memo  103/AMP  Memo  76.6 

7)  SRG  Memo  1 09/ AMP  Memo  76.7 

8)  SRG  Memo  126/ AMP  Memo  76.8 


t9.  KEY  WORDS  (Contlnuo  on  reverie  elde  II  neceeemy  and  Identify  by  block  number) 


aircraft,  equations,  kill  probabilities,  probability,  vulnerability 


20.  ABSTRACT  (Continue  on  reveree  alda  If  neceeeory  and  Identity  by  block  number) 


This  research  contribution  consists  of  a  series  of  eight  memoranda 
originally  published  by  the  Statistical  Research  Group  at  Columbia 
University  for  the  National  Defense  Research  Conriittee  in  1943  on  methods 
of  estimating  the  vulnerability  of  various  parts  of  an  aircraft  based  on 
damage  to  surviving  planes.  The  methodology  presented  continues  to  be 
valuable  in  defense  analysis  and,  therefore,  has  been  reprinted  by  the 
Center  for  Naval  Analyses  in  order  to  achieve  wider  dissemination.^ 


145 


DD  1473  COITION  OF  I  NOV  •*  IS  OSSOLITC 

S.  N  010?-  LP-  ON-  6601 


CD 

<£  0> 


<« 

§ 


O  « 


oc  f 


MEMORANDUM  FOR  DISTRIBUTION  LIST 

Sub j :  Center  for  Naval  Analyses  Research  Contribution  432 

Enel:  (1)  CRC  432,  "A  Reprint  of  'A  Method  of  Estimating 

Plane  Vulnerability  Based  on  Damage  of  Survivors' 
by  Abraham  Wald,"  July  1980 

1.  Enclosure  (1)  is  forwarded  as  a  matter  of  possible 
interest. 

2.  This  Research  Contribution  contains  a  series  of 
memoranda  written  by  Abraham  Wald  of  the  Statistical 
Research  Group  at  Columbia  University  during  World  War  II. 
Unfortunately,  this  work  was  never  published  externally, 
although  some  copies  of  his  original  memoranda  have  been 
available  and  his  methodology  has  been  employed  in  the 
analysis  of  data  from  both  the  Korean  and  Vietnam  Wars. 

It  is  published  by  CNA  not  only  as  a  matter  of  historical 
interest  but  also  because  the  methodology  is  still  relevant. 

3.  Research  Contributions  are  distributed  for  their 
potential  value  in  other  studies  and  analyses.  They  do 
not  necessarily  represent  the  opinion  of  the  Department 
of  the  Navy. 


'/2eJ  £■ 

PHIL  E.  DePOY 


PHIL  E.  DePOY 
Director 

Operations  Evaluation  Group 


jssirn  r''r 


Kris 

rrnc  ?  s 

Unp.n  •- 

Jus!.  ’  .’ica'vi.'j 


By - .  • 

DlstritJ- 

j  Ava'. 

I  * '  'Avail  ’■■■■■ 

IDist  SjecLai 


an  affiliate  of  the  University  of  Rochester 


Subj;  Center  for  Naval  Analyses  Research  Contribution  432 

DISTRIBUTION  LIST 

Department  of  the  Navy 
SNDL 

A2A  CNR 

B3  National  Defense  University  (NDU) 

B3  Armed  Forces  Staff  College  (AFSC) 

SNDL  Part  II: 

E3A  NRL 

FF38  USNA,  Nimitz  Library 

FF44  NAVWARCOL 

FKA1A  COMNAVAIRSYSCOMHQ 

FKR3C  NAVAIRTESTCEN 

FT 7 3  NAVPGSCOL 

OpNav:  Op-09BH,  Op-96,  Op-05,  Op-50W 


Other: 

Defense  Technical  Information  Center  (DTIC)  (12  Copies) 
Department  of  the  Air  Force  (SAMI) 

Air  University 

Air  Proving  Ground  Center,  Eglin  AFB 
Institute  for  Defense  Analysis  (IDA) 

Mathematics  Research  Center 
The  Mitre  Corporation 
The  Rand  Corporation 
University  of  Rochester  (10) 

Navy  Tactical’  Support  Activity 


CRC-432  /  luly  1980 


f9'ryjr^o>  / 


(y  f  A 


.  . . 

3  REPRINT  OF 

^"a'method  OF  ESTIMATING 
gLANE  VULNERABILITY  fASED  , 

— QJVQAMAGE  OF  JjURVIVORS"'T/ 

by/abraham/Wald  ;  — • 


Abraham  Wald 


(/  .////."• 

/  ‘  . - . . ; 


Operations  Evaluation  Group 


CENTER  FOR  NAVAL  ANALYSES 

2000  North  Btauragard  Stroot,  Altxandria,  Virginia  22311 


ft? ,  .oC^#Nl%4i  *  V  r  *  ;#  •  *  *  ’ 


■  Part '  iv  -  Minimum  and  Maiciaum  -; 

:;'.Vv;.v' '■  Thafe  a  PIana  will  BeCownoOt^  : 

',  Number  of  Hite  Calc ala ted  UnderSoma  Further 
¥  ■ :  on  ■■  tha . ProbabUltlaa  %#*>*  >%•• 

•’Part-V  «  Subdivision  of  the  Pi**w>  into  Several  Bqui- 

■  Area**  i^***^*  '***’  «  ••»•••.«•• . . . 


p^rt  vi  -;SonpM»9l;  '8*‘r^  *■  «  •  •  *  »  •  *  «  •  •  • 

part '  VII  ■  **  Mitotllvttout  00gKii%**  »~*  •  »•*•••  • « 

part  viii  -  Vulnerability  of  aPlana  to  Different  Types  ; 

'  of  •  Qttaa».#  *v*  »••.'♦*'••  *  ••*••  • 


‘igi1 


NffpIpjPI 


I  pl«t«  vulnerability 
lr  original  AMP  «*<£  ! 


tf'  ;  *&*•  National 

W»  ittWaahitt^ton,  0.C.  pf 


^Vvf- ■•' u;'-»V  ‘  V;^3  -¥>'j 


PART  I 


AN  EQUATION  SATISFIED  BY  THE  PROBABILITIES  THAT  A 
PLANE  WILL  BE  DOWNED  BY  i  HITS1 


INTRODUCTION 


Denote  by  P^  (i  =  1,2,...,  ad  inf.)  the  probability  that  a  plane 
will  be  downed  by  i  hits.  Denote  by  p^  the  conditional  prob¬ 
ability  that  a  plane  will  be  downed  by  the  i-th  hit  knowing  that 

the  first  i  -  1  hits  did  not  down  the  plane.  Let  Q.  =  1  -  P.  and 

ii 

qi  ~  1  "  pi  *  1,2, ... ,  ad  inf.).  It  is  clear  that 

Qi  '  qlq2***qi  (1) 

and 

Pi  "  1  “  qlq2*“qi  *  (2) 

Suppose  that  p^  and  P^  (i  =  1,2,...)  are  unknown  and  our  infor¬ 
mation  consists  only  of  the  following  data  concerning  planes 
participating  in  combats 


•  The  total  number  N  of  planes  participating  in  combat. 

•  For  any  integer  i  (i  =  0,1,2,...)  the  number  A^  of 

planes  that  received  exactly  i  hits  but  have  not  been 
downed,  i.e.,  have  returned  from  combat. 


A. 

Denote  the  ratio  jp  by  ai  (i  =  0,1,2,...)  and  let  L  be  the 
proportion  of  planes  lost.  Then  we  have 
00 

5  «i  ■  1  -  L.  (3) 

i  =  0 


*This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  85  and 
AMP  memo  76.1. 
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_  itnuk  iLfcto 


The  purpose  of  this  memorandum  is  to  draw  inferences  concerning 
the  unknown  probabilities  p^  and  P^  on  the  basis  of  the  known 

quantities  aQ,  a^,  etc. 


To  simplify  the  discussion,  we  shall  neglect  sampling  errors, 
i.e.,  we  shall  assume  that  N  is  infinity.  Furthermore,  we  shall 
assume  that 


0  <  p^  <  1  (i  =  1,2,...,  ad  inf.). 


From  equation  4  it  follows  that 


0  <  Pj;  <  1  (i  =  1,2,...,  ad  inf.). 


We  shall  assume  that  there  exists  a  non-negative  integer  n  such 
that  an  >  0  but  a^  =  0  for  i  >  n. 

We  shall  also  assume  that  there  exists  a  positive  inteqer  m  such 
that  the  probability  is  zero  that  the  number  of  hits  received  by 
a  plane  is  greater  than  or  equal  to  m.  Let  m'  be  the  smallest 
integer  with  the  property  that  the  probability  is  zero  that  the 
number  of  hits  received  by  a  plane  is  greater  than  or  equal  to 
m' .  Then  the  probability  that  the  plane  receives  exactly  m’  -  1 
hits  is  positive.  We  shall  prove  that  m’  =  n  +  1.  Since  an  >  0, 

it  is  clear  that  m'  must  be  greater  than  n.  To  show  that  rn1 
cannot  be  greater  than  n  +  1,  let  y  be  the  proportion  of  planes 

that  received  exactly  m'  -  1  hits.  Then  y  >  0  and 

y( 1  -  Pm'_i)  =  am'-l*  Since  y  >  0  and  1  -  Pm»_i  >  0,  we  have 

am, _^  >  0.  Since  a^  =  0  for  i  >  n,  we  see  that  m’  -  1  <_  n, 

i.e.,  m'<^n  +  1.  Hence,  m'  =  n  +■  1  must  hold. 

Denote  by  x^  (i  =  1,2,...)  the  ratio  of  the  number  of  planes 

downed  by  the  i-th  hit  to  the  total  number  of  planes 
participating  in  combat.  Since  m'  =  n  +  1,  we  obviously  have 

x,  **  0  for  i  >  n.  It  is  clear  that 


x .  =  L  =  1  -  a  -  a,  - 
l  o  1 


'  9 


CALCULATION  OF  x,  IN  TERMS  OF  a  ,  a ,,..., a  ,  p,,...,p 
i  o  1  n  rl  ,en 

Since  the  proportion  of  planes  that  received  at  least  one  hit  is 
equal  to  1  -  a  ,  we  have 


xa  =  Pl(l  -  aQ) 


(7) 


The  proportion  of  planes  that  received  at  least  two  hits  and  the 

first  hit  did  not  down  the  plane  is  obviously  equal  to 

1  -  a  -  a,  -  x,.  Hence, 
oil 


x  2 


-  Pj ( 1  ~ 


-  a 


"  )  • 


(0) 


In  general,  we  obtain 

x.  =  p.(l  -  aQ  -  ax  -  ...  -  a._i  -  xx  -  x2  -  ...  -  xi_i) 

(i  =  2 , 3 , . . . , n )  (9) 


i-1  ' 


Putting 

c x  =  1  -  aQ  —  ci ^  —  ...  —  a 

equation  9  can  be  written 

•••  x i _ x )  =  Pjcj  (i  =  2,3,...,n). 


(10) 


(ID 


Substituting  i-1  for  i,  we  obtain  from  equation  11 


*i-i +  + 


•  +  Xi-2)  ‘  Pi-lCi-l 


(i  =  3, 4 , . .  . ,  n ) 


(12) 


Dividing  by  pl-x,  we  obtain 


‘i-1 

’i-1 


+  (xx  +  ...  +  x^_^)  =  (i  —  3,4,...,n) 


(13) 
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•  ,  (l  -  1  \  =  qi~1  Xj  . 

1-1  \  pi-l/  pi-l  '-1 


Adding 
equation  13,  we  obtain 


to  both  sides  of 


xi  +  •**  +  xi-l  "  ci-i  ^  p: 


qi-l 


i-1 


From  equations  11  and  14,  we  obtain 


xi-l 

(i.  =  3f  4 ^  •  i n+l)  • 


(14) 


*1  +  P‘  (Ci-1  '  tel  Xi-l)  =  pi°i 


(15) 


Hence, 


piqi-l 

xi  "  pi(ci  “  ci-l)  +  ~p xi-l  {i  =  3,4,.,.,n). 


(16) 


Let 


di  *  pi(ci  ~  ci-l}  =  ”piai- 


i  i-1  (i  =  3,4, ... ,n) 


(17) 


and 


piqi-l 
1  “  Pi-i 


(i  =  3, 4, . . . ,n) .  (10) 


Then  equation  16  can  be  written  as 


*i  ”  di +  Vt-i 


{i  =  3,4, . . . ,n) .  (19) 


P2ql 


Denote  p^(l  ~  aQ)  by  d^,  “P2aj  by  d2»  and  — — —  by  t2?  then  we 
have 


X1  =  dl  and  X2  =  fc2xl  +  d2  * 


(20) 
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, .  iwr  I  MiMUteri  jUfttri'  ’ 


Prom  equations  19  and  20 ,  we  obtain 


|i  “  %  diWj+2  •“  fci+di 

EQUATION  SATISFIED  BY  <?!»••••% 

To  derive  an  equation  satisfied  by  q^  , . . . 
n 

£  x^  in  terms  of  the  quantities  t^  and 


(i  =  2, 3, . . . ,n) .  (21) 


,qn,  we  shall  express 

(i  —  1( • • • f n) • 


Substituting  i  for  i  -  1  in  equation  14,  we  obtain 


(22) 


“  ^  *" 

;i  ‘  j?i  X3  ’  °‘  ‘  Xi  ‘  C‘  '  j?!  'Vitl-'i1  +  di  • 


Hence,  in  particular 


n  I0”1 


<djtj+r"  tn>  +  dn 


(23) 


=  L  . 


n 


Since  cn  -  L  *  an,  and  since  t^  +  1  ...  =  jT"  •“  %-!'  we 


obtain  from  equation  23 


(n-1  d,  \ 

Y*  — *  q.  . .  .  q_  I  +  q„an  ,=  0  . 

h  p5  ^  7  nn-1 


(24) 


Dividing  by  q,  ...  q  and  substituting  ~P.ia.i_i  for  d.,  we  obtain 


n 


n-1 


n- 1 
-  £ 


ql  *•*  qn  ql  *•*  qn-l  j=l  Pjql  **•  qj-l 


n 


n-1 


q,  •  •  •  q  q,  • • •  q  . 

^1  ti-1 


(25) 


n-1  a.  ,  d, 

+  £  - ^ -  -  — 

j=2  qi  * ’ *  q  j-1  pl 


n 


£ 

j-i 


or 


(1  -  aQ)  =0 


(26) 


If  it  is  known  a  priori  that  qj^  =  ...  =  qn,  then  our  problem 
is  completely  solved.  The  common  value  of  q^,...,qn  is  the 
root  (between  0  and  1)  of  the  equation 


n 

h 


^4  «  1  -  a 

q] 


It  is  easy  to  see  that  there  exists  exactly  one  root  between  zero 
and  one.  We  can  certainly  assume  that  q.^  2.  q2  —  •**  —  qn*  We 

shall  investigate  the  implications  of  these  inequalities  and 
equation  26  later. 

ALTERNATIVE  DERIVATION  OF  EQUATION  26 

Let  b^  be  the  hypothetical  proportion  of  planes  that  would  have 

been  hit  exactly  i  times  if  dummy  bullets  would  have  been  used. 
Clearly  b^  >_  a^.  Denote  bi  -  a^  by  y^^  (  i  =  0, 1, 2, . . . ,n) .  Of 

n 

course,  bQ  =  aQ,  i.e.,  yQ  =  0.  We  have  £  b^  =  1.  Clearly 


-6 


yi  *  pibi  *  pi(ai  +  yl>  (i  =  l,2,...,n). 


Hence* 


1  "  ql  **•  qi 

y  x  — —  a  *  -  a .  =  — 

Qi  l  qx  ...  l  q 


ai 

-  _  a 

1  “•  qi 


Since  £  y.  =  L,  we  obtain  from  equation  28 
l-l  1 


I  ^  •  •  • 


♦  £  n  * 


1  -  a 


This  equation  is  the  same  as  equation  26.  This  is  a  simpler 
derivation  than  the  derivation  of  equation  26  given  before. 
However,  equations  21  and  22  (on  which  the  derivation  of  equation 
26  was  based)  will  be  needed  later  for  other  purposes. 


As  mentioned  before,  equation  29  leads  to  a  solution  of  our 
problem  ii  it  is  known  that  =  ...  «  qn»  In  the  next 

memorandum  (part  II)  we  shall  investigate  the  implications  of 
equation  29  under  the  condition  that  q^  >,  q.>  >_  ...  _>  qn« 


NUMERICAL  EXAMPLES 


N  is  the  number  of  planes  participating  in  combat.  Aq,  A^,  a.,. 


...,An  are  the  number  returning  with  no  hits,  one  nit,  two  hits. 


...,n  hits,  respectively.  Then 


_  Ai 
di  N 


(i  =  0,1, 2, ...,n) 


i.e.,  is  the  proportion  of  planes  returning  with  i  hits.  Ti 


x 

computations  below  were  performed  under  the  following  two 
assumptions : 


mmtp 


•  The  bombing  mission  is  representative  so  that  there  is  no 
sampling  error. 

•  The  probability  that  a  plane  will  be  shot  down  does  not 
depend  on  the  number  of  previous  non-destructive  hits. 


N 

S 

400 

Ao 

= 

320 

then  a^ 
o 

= 

o 

CO 

• 

A1 

= 

32 

al 

= 

.08 

A2 

= 

20 

a2 

= 

.05 

A3 

= 

4 

a3 

= 

.01 

A4 

a 

2 

a4 

= 

.005 

A5 

s 

2 

a5 

= 

.005 

We  assume  q^  =  q^  ■  . . .  »  q5  *  q^,  where  q^  is  the  probability  of 

a  plane  surviving  the  i-th  hit,  knowing  that  the  first  i  -  1  hits 
did  not  down  the  plane. 

Then  equation  26, 


reduces  to 
n  a . 


—  =  1  "  ao  ' 
J 


S.  3  =  1  **  ao  * 
i=i  q 


Substituting  values  of  a. 

l 


.08  .05  .01  .005  .005  _n 

+  2  +  3  +  4+  5  ~  • 

q  q  q  q  q 


WPS*"1  »'»  ■igi.«>-<g"y  -.--  -■»  •-  '&***■  ***■-  >•<*">-* 


or 


.  200q5  -  .080q4  -  .  050q3  -  .OlGq2  -  .005q  -  .005  =  0. 


The  Birge-Vieta  method  of  finding  roots  described  in  Marchant 
Method  No.  225  is  used  to  solve  this  equation  (table  1).  We  find 
q  =  q^  -  .851,  p^  =  .149  where  is  the  probability  of  a  plane 

being  downed  by  the  i-th  hit,  knowing  that  the  first  i  -  1  hits 
did  not  down  the  plane. 

equals  the  ratio  of  the  number  of  planes  downed  by  the  i-th 

hit  to  the  total  number  of  planes  participating  in  combat.  Using 
equation  9 


$ 

1 

xi  =  pi<x 

-  a 

o 

U1 

•  • 

•  ™  a^  ^  *  x  2  "*  x  ^  —  •  •  •  * 

xi-l) 

& 

t 

t 

(i  =  2,3,. 

• .  ,n) 

i 

£ 

for  n  =  5,  we  obtain 

!' 

X1  =  PXU  -  ao) 

“  • 

030 

8 

• 

x2  =  P2(l  -  aQ 

-  al 

-  xx) 

013 

i 

x3  =  p3(l  -  aQ 

~  al 

-  a2 

_  X1 

-  x2)  =  .004 

l 

1 

X4  =  P4(1  “  ao 

-  al 

a2 

d3 

-  xi  -  x2  -  x3)  =  .002 

?. 

y 

x 5  »  -  a0 

al 

a2 

d3 

-  a4  -  x3  -  x2-  x3  -  x4)  = 

.001 

example  2 :  Let 

a 

o 

=  •  3 1 

al  = 

•2,  a2  =  .1,  a3  =  .1,  a^  = 

.05, 

a_  =  .05.  Then 
5 

the 

following 

results  are  obtained:  q  = 

.  07, 

P=  1  -  q  =  .13 

'  X1 

=  .09 

,  x2 

=  .05,  x3  =  .03,  x4  =  .02, 

and 

x  5  =  .01. 


The  value  of  q  in  the  second  example  is  nearly  equal  to  the  value 
in  the  first  example  in  spite  ot  the  fact  that  the  values  a^^ 

(i  =  0,1,..., 5)  differ  considerably.  The  difference  in  the 
values  a^  in  these  two  examples  is  mainly  due  to  the  fact  that 

the  probability  that  a  plane  will  receive  a  hit  is  much  smaller 
in  the  first  example  than  in  the  second  example.  The  probability 
that  a  plane  will  receive  a  hit  has,  of  course,  no  relation  to 
the  probability  that  a  plane  will  be  downed  if  it  receives  a 
hit . 
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TAULE  1 


ftMUM  q 

-  «  -  y, 

.200 

-.000 

-.050 

-.  0  10 

-.005 

-.005 

♦  .200 

♦  .  120 

♦  .070 

♦  .060 

♦  .055 

.200 

♦  .  120 

♦  .070 

♦  .060 

♦  .055 

♦  .050 

♦  .200 

♦  .320 

♦  .  390 

♦  .450 

.200 

♦  .  320 

♦  .390 

♦  .450 

♦.505  - 

*1 

A 

y  _  -  y .  -  ~  »  1  -  .0990  -  .9010 

2  1  A 


Assume  q 

-  .9010  - 

*2 

•  2000 

-.0800 

-.0500 

-.0100 

-.0050 

-.0050 

♦  .  1802 

♦.0903 

♦.0363 

♦.0237 

♦•0168 

.2000 

♦  .  1002 

♦.0403 

♦.0263 

♦  .0  187 

♦.0118  «  B 

♦ . 1802 

♦.2526 

♦.2639 

♦.2615 

o 

.2000 

♦.2904 

♦. 2929 

♦. 2902 

♦.2002  - 

", 

B 

y3  -  y2  - 

--  -  .9010 

1 

-  .042113 

-  .858887 

Assume  q 

-  .858887 

.200000 

-.080000 

-.050000 

-. 0  10000 

-.005000 

-.005000 

♦. 171777 

♦.070826 

♦.024758 

♦.012675 

♦.006592 

.200000 

♦. 091777 

♦.028826 

♦  .  0  14758 

♦.007675 

♦.001592  -  C 

♦  .  171777 

♦.226363 

♦.219179 

♦.200925 

o 

. 200000 

♦.263554 

♦  .255  189 

♦.233937 

♦.200600 

*  C1 

C 

y 

4  “  y3  “  C 

-  -  .858887 
1 

-  .007632 

-  .851255 

Assume  q 

-  .851255 

*  y4 

.2000000 

♦ .000000 

-.050000 

-.0  10000 

-.005000 

-.005000 

♦. 17025 1 

♦.076027 

♦.022837 

♦.010928 

♦.005046 

.2000000 

♦.090251 

♦.026627 

♦.012037 

♦. 005920 

♦.000046  -  D 

♦ . 170251 

♦.22 1754 

♦  .2  1  1606 

♦. 191058 

o 

•2000000 

♦.260502 

♦•24853 1 

♦.224443 

♦. 196986 

-  D. 

y5  •  r4  -  J~  ■  -BSUSi  -  .0002)4  -  .85  )02  1 
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PART  II 


MAXIMUM  VALUE  OF  THE  PROBABILITY  THAT  A  PLANE  WILL  BE  DOWNED 

BY  A  GIVEN  NUMBER  OF  HITS1 

The  symbols  defined  and  the  results  obtained  in  part  I  will  be 
used  here  without  further  explanation.  The  purpose  of  this 

i 

memorandum  is  to  derive  the  least  upper  bound  of  X .  =  £  x.  and 

j  =  1 

that  of  P^^  (i  =  l,...,n)  under  the  restriction  that 
qi  >_  q 2  2.  •••••  2.c^n* 

First,  wo  shall  show  that  Xi  is  a  strictly  increasing  function 
of  Pj  for  j  £  i  •  Let  us  replace  p^  by  p^  +  A  (A  >  0)  and  let 
us  study  the  effect  of  this  change  on  x^,...,x^.  Denote  the 
changes  in  x^,...,x^  by  A^,...,A^,  respectively.  Clearly, 

A^  =  ...  =  A^  =  0.  It  follows  easily  from  equation  9  that 

Aj  >  0  and 

Aj  +  1  =  "V1  Aj  * 

Hence, 

h  +  4j+i  '  (1  -  p:*i’  h  > 

Similarly,  we  obtain  from  equation  9 

A j  +2  =  -pj+2(Aj  +  Aj  +  1}  =  “pj+2(1  “  pj  +  lJ  Aj  * 

Hence , 

4j  +  V*  *  4jt2  ■  (1  -  Pj+2>  11  -  Pj+1>  4j  >  °' 

In  general 

4j  +  4j+1  ♦  ...  +  Aj+k  -  U  -  Pj+1>  ...  (1  -  Pj+k>  Sj  >  0 

(k  =  1, • • • , i-j ) 

Hence,  we  have  proved  that  X^  is  a  strictly  increasing  function 
of  Pj  ( j  *  1, . . . , i) . 

1This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  E7  and 
AMP  memo  76.2. 
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On  the  basis  of  the  inequalities  we  shall  derive  the 

least  upper  bound  of  X^ .  For  the  purpose  of  this  derivation  we 

shall  admit  0  and  1  as  possible  values  of  (i  =  l,...,n),  thus 

makinq  the  domain  of  all  possible  points  (p^,...,pn)  to  be  a 

closed  and  bounded  subset  of  the  n-dimensional  Cartesian  space. 
Since  X^  is  a  continuous  function  of  the  probabilities  p^,  p^, 

etc.  (X  is  a  polynomial  in  p^,...,p^),  the  maximum  of  X^  exists 

and  coincides,  of  course,  with  the  least  upper  bound.  Hence, 
our  problem  is  to  determine  the  maximum  of  X^. 

First,  we  show  that  the  value  of  X.  is  below  the  maximum  if 

i 

p  >  (>..  Assume  that  p  >  p.  and  let  k  be  the  smallest  positive 
rn  1 l  'n  1  l 

inteyer  for  which  p.>  p .  .  Obviously  k  >  i.  Let  p'  =  p  (1  +  e) 

K  1  j  j 

tor  j  =  l,...,k-l,  and  p!  «  p^  (1  -  n)  for  3  =  k,k+l,...,n, 

where  e  >  0  and  n  is  a  function  n (  e )  of  e  determined  so  that 
n 

£  x'  =  L  (x.j  is  the  proportion  of  planes  that  would  have  been 

j"l  3  3 

brought  down  with  the  j-th  hit  if  pj,...,p^  wore  the  true 

probabilities).  Since  Xr  (r  =  l,...,n)  is  a  strictly  monotonic 

lunction  of  p^,...,p  ,  it  is  clear  that  tor  sufficiently  small 

such  a  function  h(  e  )  exists.  It  is  also  clear  that  tor  suffi¬ 
ciently  small  e  the  condition  p'.<  p'  <  ...  <  p'  is  fulfilled. 

Since  pj  >  pj  (3  =  l,...,i),  we  see  that  X|  >  X^  (X^  does  not 

depend  on  p^.  for  r  >  i).  Hence,  we  have  proved  that  if 

Pj,...,Pn  is  a  point  at  which  becomes  a  maximum,  we  must  have 

Pi  =  Pi  +  JL  =  ...  =  P„. 

How  we  shall  show  that  if  X^  is  a  maximum  then  p^  =  p(  =  ...  =  p^ 
For  this  purpose  assume  that  p^  >  p^  and  we  shall  derive  a  con¬ 
tradiction.  Let  j  be  the  greatest  integer  for  which  p^=  p^. 
Since  p^>  p^,  we  must  have  3  <  i.  Let  p^.  =  pf(l  +  £  )  for 

r  *  l,...,j  and  p£  =  pr(l  -  n  )  for  r  =  j+l,...,i,  where 

i  i 

e  >0  and  n  is  determined  so  that  Y  x^  =  Y  x.  .  Then  for  the 

k^l  k  k= 1  k 

probabilities  pj,...,p[,  p^+^,...,pn  the  proportion  of  lost 


planes  is  not  changed,  i.e.,  it  is  equal  to  L.  Now  let  ■  p!^ 

for  r  >  i.  Then  the  proportion  L*  of  lost  planes  corresponding 
to  pj,...,p^  is  less  than  L.  Hence,  there  exists  a  positive 

A  so  that  the  proportion  L"  of  lost  planes  corresponding  to  the 
probabilities  P£  =  (1  +  A )  is  equal  to  L.  But,  since  p£  >  p^. 

i  i  i 

(r  ■  l,...,i)  we  must  have  Y]  x?  >  V  x\  =  V  x..  Hence,  we 

j=l  3  j=l  3  j=l  3 

arrived  at  a  contradiction  and  our  statement  that  p.  =  Pj  =  . . .  “ 

p^  is  proved.  Thus,  we  see  that  the  maximum  of  is  reached 

when  p,  *  p.  *  ...  *  p  . 

1  2  n 

LEAST  UPPER  BOUND  OF  Pj 

Now  we  shall  calculate  the  least  upper  bound  of  P^ .  Admitting 
the  values  0  and  1  for  p^ ,  the  maximum  of  P^^  exists  and  is  equal 
to  the  least  upper  bound  of  P^ .  Since  Pi  =  1  -  q^  ...  q^ 
maximizing  P^  is  the  same  as  minimizing  q^  ...  q^.  We  know  that 
ql'*.*,qn  are  sub3ect  to  the  restriction 


IP 


j-1  ql 


a  1  -  a. 


(30) 


Let  q°,, 


. . ,qn  be  a  set  of  values  of  q^,...,qn  (satisfying 


equation  30)  for  which  q^  . . .  qj  becomes  a  minimum.  First,  we 


show  that 


o 

qi+l 


‘n 


Suppose  that  q°  <  q° 


Consider  the  set  of  probabilities  q£  *  q°  for  r  <  1  and  q£.  *  q? 


for  r  >  i .  Then 


3  =  1 


q r 


fi 


<  1  -  a. 


-13- 


Hence,  there  exists  a  positive  factor  A  <  1  so  that 


where  q?  =  Aq|  (i  * 


Then 


qIq2 


qi  <  q!q5 


in  contradiction  to  our  assumption  that  q°  ...  cj?  is  a  minimum. 

Hence,  we  have  proved  that  q?  =  ...  =  q°  . 

1  n 

Now  we  show  that  there  exists  at  most  one  value  j  such  that 
1  >  q°  >  q?  .  Suppose  there  are  two  integers  j  and  k  such  that 

1  >  q?  >  q°  >  q?  .  Let  j'  be  the  smallest  integer  for  which 
j  —  x  1 

q?,  =  q?  and  let  k'  be  the  largest  integer  for  which  q°  *  q°  . 
J  J  K  K 

Let  qjt  =  {1  +  e)  q°,,  qk,=  q£,  (£  >0),  and  qr  =  q° 

for  r  fi  j '  ,  /  k'  .  Then 

„  „  n  a 

ql  *•*  qi  *  ql  ***  qi  and  2*  Z - ~  <  1  “  a0  • 

r=l  ^1  *r 

Hence,  there  exists  a  positive  factor  A  <  1  such  that 


where  q*  =  Xq^  .  but  q£  ...  q*  <  q^  ...  qi  =  q°  ...  q°,  which 

contradicts  the  assumption  that  q°  ...  q°  is  a  minimum.  This 
proves  our  statement. 


It  follows  from  our  results  that  the  minimum  of  q^  is  the  root  of 
the  equation 

n  a 


E  -f  -  *  - 

r=l  q 


a. 


(32) 


Now  we  shall  calculate  the  minimum  of  2 •  First,  we  know  that 
q^  =  q2  (i  >_  2)  if  q ^q 2  be  a  minimum.  Hence,  we  have  to  minimize 
q^q^  under  the  restriction 


1  L  +  2  +  3  +  .  n  \ 

q"T  (a  1  +  “  +  ~2  +  +  “n^I  ~ 

1\  q2  q2  q2  / 


=  1  -  a 


(33) 


Using  the  Lagrange  multiplier  method  we  obtain  the  equations 


<3  ;  " 


5  fi +  7  *7  +  •••  0 

ql  \  q2  q2  q2  / 


(34) 


(Lagrange  multiplier  =  A) 


±-(ll  +  111  +  +  (n  '  1>an\  . 

“q,(  2  3  +  *  ’  *  n  1“ 

1  \  q2  q2  q2  / 


(35) 


Because  of  equation  33,  we  can  write  equation  34  as  follows: 
'  qlq2 


-t,  u 


Substituting  for  A  in  equation  35,  we  obtain 


1  /a 2  .  2a3  .  Ja4  .  .  <n  ’  1)an\  „ 

ql  "  T~-  a  I  ~  +  ~T  +  "T-  - n^T -  )  =  0 

\q2  q2  q2  q2  / 


(36) 
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or 


1  /a2  .  2a3  (n  -  1,an^ 

qi  1  -  ao  I  2  **•  +  n-1  /  * 

\q2  q2  q2  / 


On  the  other  hand,  from  equation  33  we  obtain 


qi  "  T 


L  +  ^  +  -U  + ...  +  . 

°\  q2  q2  q2  / 


Equating  the  right-hand  sides  of  equations  37  and  38,  we  obtain 


(n  -  2)a 


“  a ^  ■  0. 


It  is  clear  that  equation  39  has  exactly  one  positive  root.  The 
root  is  less  than  or  equal  to  1  if  and  only  if 


a^  +  2a^  +  3a5  +  ...  +  (n  -  2)an  <  a^  . 


Equations  38  and  39  have  exactly  one  positive  root  in  q^  and  q^. 

We  shall  show  that  if  the  roots  satisfy  the  inequalities  1  >  >  q^, 

then  for  these  roots  q^q^  becomes  a  minimum.  We  can  assume 
that  2  <  n,  since  the  derivation  of  the  minimum  value  of  ...  qn 
will  be  given  later  in  this  memorandum.  It  is  clear  that  for  any 

al 

value  q^  >  ^  -  equation  38  has  exactly  one  positive  root  in 

o 

q2«  Denote  this  root  by  Hence,  Mq^)  is  defined  for 


00 


lim  <f>  (q^ )  =  + 


Hence  (assuming  a^  >  0) 


lim  ip  (q^ )  =  +  oo 


’l  *  1=31 


"1 


where  <1*^)  =  q1  (J^)  . 

It  is  clear  that  lim  $(q.)  =  0.  Since  a  >  0,  it  follows  from 

q^oo  n 

equation  38  that  q^ (  0 ( q^ ) ] °  1  has  a  positive  lower  bound  when 
q,+  00  .  But  then,  since  n  >  2,  lim  q,  0  (q,  )  =  +  «.  From 

the  relations  lim  (p(q..)  =  lim  ^(q.  )  *  +  »  it  follows 

O 

that  the  absolute  minimum  value  of  ^(q^)  is  reached  for  some 
positive  value  q^  .  Since  equations  38  and  39  have  exactly  one 
positive  root  in  q^  and  q^,  the  absolute  minimum  value  of  0 ( q ^ ) 

must  be  reached  for  this  root.  This  proves  our  statement  that 
if  the  roots  of  equations  38  and  39  satisfy  the  inequalities 
1  ^  q^  ^  q2,  then  for  these  roots  q^q2  becomes  a  minimum  con¬ 
sistent  with  our  restrictions  on  q^  and  q2»  If  1  q^  q2  is 
not  satisfied  by  the  roots  of  equations  38  and  39,  then  q^  is 
equal  either  to  1  or  to  q2  and  the  minimum  value  of  q^q2  is 

2 

either  $(1)  or  q  ,  where  q  is  the  root  of  the  equation 
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n  ar 

T  —  =  1  -  a 
r=l  qr  ° 


Mow  we  shall  determine  the  minimum  of  ...  q^(2<i<n). 
First,  we  determine  the  minimum  M.^  of  ...  q^  under  the  re¬ 
striction  that  =  q^.  Thus,  we  have  to  minimize  q^q2  ^  under 
the  restriction  that 


al  2  3 

— i  4.  t —  4.  — 

ql  qlq2  qlq2 


+  ...  + 


qlq2 


=  1  -  ao. 


(40a) 


Using  the  Lagrange  multiplier  method,  we  obtain 


*l-  tk . 


n  \  i-1  X  ,  , 

q  qn_1  )  =  q'2  ~  qi  ”  a°  "  ‘ 5 

qlq2  /  1  (41) 


<i  -  llq^'*2  -  ^  (V  + 


(n  -  1 ) a 


+  .  . .  + 


=  0  . 


\q2  q2  q2  /  (41a) 

qlq2_1 

Substituting  y— — —  for  X  (the  value  of  X  obtained  from 


equation  41),  we  obtain 


u  -  uq,  -rh 


a  2a 

—  +  —  +  •••  + 


(n  -  l)ar 


°  \q2  q2 


=  0  . 


From  equation  40a 


(1  -  1)qi  '  K-r 


ai  + 


=  0  .  (43) 


Httf'ilWTi  ifr**~f* 
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Prom  equations  42  and  43,  we  obtain 


(i 


(i  -  2)a, 


l)a, 


(i  -  n)a 


n 


_n-l 


=  0. 

(44) 


From  Descartes'  sign  rule  it  follows  that  equation  44  has  exactly 
one  positive  root. 


Let  q^  =  q°  and  -  q°  be  the  roots  of  the  equations  43  and 

44.  If  1  >  q°  >  q°  »  then  M.^  =  q°(q°)1_1  .  if  i  >.  q°  >  q°  does 

not  hold,  then  is  either  (q * ) 1  or  (q")1-1,  where  q'  is  the 

root  of  the  equation 


n 

£ 


=  1  -  a 


]=1  (q')3 

and  q"  is  the  root  of  the  equation 


al  + 


n 


( q " ) 2 


<q") 


n  »n-l 


=  1  -  a. 


(45) 


(46) 


Let  (r  =  2,...,i-l)  be  the  minimum  of  q^  ...  q^^  under  the 

restriction  that  q.  -  ...  »  q  ,  =  1  and  q  , ,  =  q .  .  Then  M.  can 

1  r-1  r+1  l  lr 

be  calculated  in  the  same  way  as  M^;  we  have  merely  to  make  the 
substitutions 


* 


n 

S 

n  -  r  + 

1 

* 

ao 

* 

3 

ao  +  al 

+  •  •  • 

+  ar-l 

* 

a3 

* 

“ 

a j+r-1 

( j  =  1, . . . ,n  ) 

* 

qi 

qj+r-l 

( j  =  lr . • . ,n  ) 

* 

l 

* 

i  -  r  + 

1, 

*  * 

and  we  have  to  calculate  the  minimum  of  ...  q^*  .  Thus,  we 

have  to  solve  the  equations  corresponding  to  equations  43  and  44, 
i.e.,  the  equations 


if  -  i 


19- 


(43*) 


and 

*  *  *  *  *  *  * 

*  *  (i  -  2)a„  (i  -  3)a-,  (i  -  n  )a  * 

(i  -  l)a  +  - -  +  - 5-i  +  ...  +  - —  =  0. 

X  *  *  *  *  n*— 1 

q2  (q2>  (q2>  (44*} 

*  * 

Let  q^  =  v^  and  q2  =  be  the  positive  roots  of  the  equations 

* 

43*  and  44*.  If  1  ^  v^  _>  v^t  then  M^r  =  v^v*  ^  ^  2.  1  v2 

* 

does  not  hold,  then  Mir  is  equal  to  either  ( v • ) 1  or 
* 

(v")1  -1,  where  v'  is  the  positive  root  of  the  equation 

*  * 


and  v"  is  the  positive  root  of  the  equation 


* 

*  a2 
a.  +  - 

* 

♦  2 

+  ...  +  1 

* 

a  * 
n* 

* 

=  1  -  a 

(46*) 

1  v" 

(v-)2 

(V" 

j  n*-l 

o 

The  minimum  t 

h  of  ql 

...  qt  (i  = 

2,3,. 

..,n-l)  is  equal 

to  the 

smallest  of  the  i  -  1  values  M. ,,..., M.  .  .  . 

If  1“ X 

Now  we  shall  determine  the  minimum  of  q^  ...  qn  .  We  show  that 

the  minimum  is  reached  when  q^  =  ...  =  qn_^  =  1.  Suppose  that 

this  is  not  true  and  we  shall  derive  a  contradiction^  Let  j  be 
the  smallest  integer  for  which  q^  <  1  (j  <  n) .  Let  q_.  =  (1  +  e)cj 

(e  >  0),  qn  *  —  -  ,  and  qf  =  qf  for  all  r  /  j,  /  n. 

Then  q, ...  q  ...  *  q,...  q  and 
1  n  1  n 
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E  — 

r=l  qx 


<  1  -  a. 


Hence,  there  exists  a  positive  X  <  1  such  that 


E  ■«— — *  -  1  -  V 

r=l  q1  ...  qr 


where 


*  — 


But  then  q,...  q^<q,...q^*q,  ...  q  in  contradiction  to 
1  n  1  n 

the  assumption  that  q^  ...  qn  is  a  minimum.  Hence,  we  must  have 
q^  ■  ...  =  qn-1  =  1.  Then,  from  equation  26  it  follows  that  the 
minimum  value  of  q^  ...  qn  is  given  by 


1  *  a  -  a .  “  ...  —  a  . 
o  1  n— 1 


If  i  >  1  but  <  n,  the  computation  of  the  minimum  value  of  q^  ...  c 

is  involved, . since  a  large  number  of  algebraic  equations  have  to 
be  solved.  In  the  next  part  we  shall  discuss  some  approximation 
methods  by  means  of  which  the  amount  of  computational  work  can  be 
considerably  reduced. 


PART  III 


APPROXIMATE  DETERMINATION  OF  THE  MAXIMUM  VALUE  OF  THE  PROBABILITY 
THAT  A  PLANE  WILL  BE  DOWNED  BY  A  GIVEN  NUMBER  OF  HITS* 

The  symbols  defined  in  parts  I  and  II  will  be  used  here  without 

further  explanations.  We  have  seen  in  part  II  that  the  exact 

determination  of  the  maximum  value  of  P.  (i  <  n)  involves  a  con- 

1 

siderable  amount  of  computational  work,  since  a  large  number  of 
algebraic  equations  have  to  be  solved.  The  purpose  of  this 
memorandum  is  to  derive  some  approximations  to  the  maximum  of  P^ 

which  can  be  computed  much  more  easily  than  the  exact  values. 


Let  us  denote  the  maximum  of  P^  by  P°  and  let  0°  =  1  -  P°  • 

Thus,  Q°  is  the  minimum  value  of  .  Before  we  derive  approxi¬ 
mate  values  of  P°  (or  Q°)  we  shall  discuss  some  simplifications 
that  can  be  made  in  calculating  the  exact  value  P°  (or  Q°) 

assuming  1  <  i  <  n.  We  have  seen  in  part  II  that  Q?  is  equal  to 
the  smallest  of  the  i  -  1  values  M. .,..., M,  .  We  shall 

make  some  simplifications  in  calculating  M^  (r  =  l,...,i-l). 


For  this  purpose  consider  the  equation 

a_ 

-  +  ...  +  — 

U  UV  UV 


a  a  , . 


"  X  3  *  3  *  “  •••  “cl  i  ■ 

n-r  o  1  r-1 


(47) 


It  is  clear  that  for  any  value  u  >  ■_ 


ao  "  ***  ar-l 


,  equation 


47  has  exactly  one  positive  root  in  v.  Denote  this  root  by  <f>r(u) 


Thus,  ♦  (u)  is  defined  for  all  values  u  > 


1  "  ao“  •••  ar-l 


In  all  that  follows  we  shall  assume  that  a^  >  0  (i  =  l,...,n). 
We  shall  prove  that 


U“  (u  kr(u)l  i“r)  " 

\  l  J  / 


=  +  oo 


(4U) 


J-This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  Uti  and 
AMP  memo  76.3. 
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r 


L 


and 


lim 

U  -+oo 


S  +00 


(49) 


It  follows  easily  from  equation  47  that  if  u  -*•  ^ 


a 

o 


-  a 


r-1 


then  ♦r(u) 


-*•  + 1 


Since  i  >  r,  we  see  that  equation  48  must 


hold.  It  follows  easily  from  equation  47  that  lim  4  (u)  =  0. 

U=+  oo  r 


Vie  also  see  from  equation  47  that  if  u>«,  the  product 

must  have  a  positive  lower  bound.  Equation  49  follows 
and  the  fact  that  lim  $r(u)  -  0. 


u[*r(uj)n  C 

from  this 


U  oo 


We  have  seen  in  part  II  that  equations  43*  and  44*  have  exactly 


one  positive  root  in  the  unknowns,  q^  and  q2*  Let  the  root  in 


*  .  o 

q,  be  u .  . 
^1  xr 


Then  the  root  in  q_  is  equal  to  4  (u.  ). 

2  r  xr 


From  equations  48  and  49  it  follows  that  u  ^4r(u)J  i-r  i 


xs 


strictly  decreasing  in  the  interval  j— — 


ao  “  “  ar-l 


<  u  <  u-r. 


and  is  strictly  increasing  in  the  interval  u°r  <  u  <  +  < 


Denote  by  u£  the  positive  root  of  the  equation 


ar  **r+l 

—  +  j-  + 

u  u 


n-r+1 


=  1  -  a_  - 


-  a 


r-1* 


(50) 


It  is  clear  that  u£  <  1  and  4r(up  =  u£ 


The  value  Mir  is 


equal  to  the  smallest  of  the  three  values 


ur  [Vur']1"r'  [♦r(1)]1'r-  and  uir  [♦r(u!r,]*'C' 


A  simplification  in  the  calculation  of  can  be  achieved  by 


the  fact  that  in  some  areas  M^r  can  be  determined  without 


calculating  the  value  u^  .  We  consider  three  cases. 


Case  At 


“i  [♦r<uP]i"r  <  [♦r«1»]i'r 
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In  this  case. 


Mir  =  ur  [<J>r(ur,J1~r  if  du  u  [*r(u)]  ^  >  0  for  u  =  u£ 


and 


M  •  =  u  . 

lr  ir 


Case  B: 


In  this  case. 


[♦r<uir']1*r  i£  dU  u[*r(uj1"r  <  0  for  u  *  ur' 


Mir  =  4>r(  D  i"C  if  u[Vu^  <  0  for  u  =  1 


and 


M.  =  u° 
ir  xr 


Case  C:  u' 


In  this  case. 


M .  =  u . 

ir  xr 


[Vulrl]i"r  i£  3^  u[*r(u>]1‘r  >  0  for  u  -  1. 

[<’'ur>]i'C-  [♦»>]  i‘r  • 


-r 


We  can  easily  calculate  the  value  of 
and  u  =  1.  In  fact,  we  have 


du 


u  [v>] 1_r 


for  u  =  u ' 
r 


la  j*r<urj1~r  *  [1>r(u)]i’r  +  ^  -  r>u  ^r(u>]1 


,  d$>  ( u ) 
i-r-1  rr 


"3u" 


(51) 


d*  (u)  dy 

and  — ga —  =  can  be  obtained  from  equation  47  as  follows, 
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a  a  i  a 

Denote  — —  +  — - —  +  ...  +  — n_r  by  G(u,v).  Then 

u  uv  uv 


d<|>r(u) 

"“du 


dv 

du 


&  G|u'v) 
G(u,v) 


(52) 


1 

u 


r+1 


uv 


n 


uv 


n-r 


(n  -  r)a 


1 

u 


r+1 


Ta 


r+2 


n) 


n-r+1 


-  "  ao  "  dl 


-  a  ) 
r-1 


r+1 


2a 


r+2 


(n  -  r)af 
n-r+1 


^v  =  *r(u)^  . 


on  the  uasis  of  equations  51  and  52»  we  can  easily  obtain  the 


value  of 


d  L  -  .  i-r 
cTuu  rr(u 


for  u  =  u^.  and  u  ■  1  if  u£.  and  ♦  r ( 1 ) 


have  been  calculated.  If  u  a  u^  ,  then  $r(u)  =  v  =  u^.;  if 
u  =  1,  then  v  =  ♦r(l). 

Since  ♦f(l)  is  equal  to  the  root  of  the  equation  in  v 


a  , ,  a 

,  r+1  .  .  n 

r  v  v 


=  1  "  d .  •  3-1  *  *  *  •  ^  3  i  § 

n-r  o  1  r-i 


it  follows  from  equation  50  that 


♦  r(l) 


r+1  * 


(53) 


\ 


* 
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Thus,  for  carrying  out  the  investigations  of  cases  A,  B,  and  C 
for  r  =  l,...,i-l,  we  merely  have  to  calculate  uj,...,up 

If  we  want  to  calculate  Q°  for  all  values  i  <  n,  then  it  seems 
best  to  compute  first  the  n  quantities  uj,...,u^. 

Since  u£  =  $r(up  and  <$>r(l)  =  u£.+  i,  we  can  say  that  is 

the  smallest  of  the  three  values 

Kr+1'  M*".  -  uir  [V“T, SY-'  • 

Since  Q?  is  equal  to  the  minimum  of  the  i  -  1  values, 

M. M.  .  .,  we  see  that 

X  JL  1  f  1“1 

Q°  <  t.,  (54) 

where 

(u*)1"1,...,  (u[_x)2,  u?J  .  (55) 

be  expected  that  Q°  will  be  nearly  equal  to 

t^.  Thus,  t^  can  be  used  as  an  approximation  to  Q?  .  In  order 
to  see  how  good  this  approximation  is,  we  shall  derive  a  lower 
bound  z^  for  Q°  .  If  the  difference  t^  -  z^  is  small,  we  are 

certain  to  have  a  satisfactory  approximation  to  Q°  .  If  t.  -  z^ 

is  large,  then  t^  still  may  be  a  good  approximation  to  0°,  since 

it  may  be  that  z^  is  considerably  below  Q°  . 

To  obtain  a  lower  bound  z^  of  Q°,  denote  by  y^  (j  =  0,1,..., i-1) 

the  proportion  of  planes  (number  of  planes  divided  by  the  total 
number  of  planes  participating  in  combat)  that  would  be  downed 
out  of  the  returning  planes  with  j  hits  if  they  were  subject  to 
i  -  j  additional  hits.  Then 

Pi  "  yo  +  yl  +  *•*  +  yi-l  +  xi  +  x2  +  •••  +  xi  •  (56) 


ti  =  Min  (up1, 


If  n  is  large,  it  can 
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i 


It  is  clear  that  ajPi>  yj  (j  *  0,1,..., i-1)  and  consequently 
(ao  +  ai  +  ...  +  a.^)  P.  >  yQ  +  yx  +  ...  +  y.^  . 


Hence , 


yo  *  yl  *  —  +  yi-l  .  p 

a^  +  a ,  +  . . .  +  a  .  ,  1 

o  1  1-1 


Equation  56  can  be  written 


pi  =  (ao  +  •••  +  *!_!> 


yo  *  yl  *  4  yi-l 

ao  +  *•*  +  ai-l 


X.  +  ...  +  X  . 

+  (1  -  aQ  -  ...  -  I  _  a  _  >t  _  _'a  • 

O  1  —  1 

yO  +  •••  +  ^i-l 

Hence,  P.  is  a  weighted  average  of  — — 7 - r— - -  and 

1  aQ  +  . . .  +  ai-JL 


xi  +  ...  +  xA 
-  aQ  -  ...  -  a._1 


Then,  from  equation  57  it  follows  that 


x^  +  ...  +  x^ 

Pi  <  1  -  ao  -  aA  -  ...  -  ai_A 


Since  y^  >  0,  we  obtain  from  equations  56  and  59 


xA+...+xA<P^< 


x^  +  . . .  +  x^ 

-  ao  -  ...  -  ai_1 


Hence , 


x.  +  ...  +  x. 

1  “  1  -  a-  ...  -  a  <  Qi  <  1  “  {X1  +  *•*  +  xi)  * 
o  1—1 
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In  part  II  we  have  calculated  the  maximum  value  of  x,  +  ...  +  x.. 

1  i 

Denote  this  maximum  value  by  .  Then  a  lower  bound  of  Q°  is 
given  by 


NUMERICAL  EXAMPLE 

The  same  notation  will  be  used  as  in  the  numerical  examples  for 
part  I.  q^  is  the  probability  of  a  plane  surviving  the  i-th  hit, 

knowing  that  the  first  i  -  1  hits  did  not  down  the  plane.  Then 
the  probability  that  a  plane  will  survive  i  hits  is  given  by 

Qi  "  ^2  qi  ‘ 

In  part  I  it  was  assumed  that 

qi  =q2  =  *..  =qi  =qo  (Say) ' 

which  is  equivalent  to  the  assumption  that  the  probability  that  a 
plane  will  be  shot  down  does  not  depend  on  the  number  of  previous 
non-destructive  hits.  Under  this  assumption 


The  example  below  is  based  on  the  assumption  that 
ql  >  q2  >  •••  >  <V 

i.e.,  the  probability  of  surviving  the  i  +  1  hit  is  less  than  or 
equal  to  the  probability  of  surviving  the  i-th  hit.  In  this 
case,  it  is  not  possible  to  find  an  explicit  formula  for  ,  but 

a  lower  bound  can  be  obtained.  That  is,  a  value  of  Q.  can  be 

found  such  that  the  actual  value  of  must  lie  above  it.  The 

o 

greatest  lower  bound  is  denoted  by  .  Hence ,  we  have 

Q-  <  Qi  • 
if 
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P?  is  the  least  upper  bound  o£  ;  that  is,  the  probability  of 

being  downed  by  i  bullets  cannot  be  greater  than  P?. 

Since  the  computation  of  the  exact  value  of  Q°  is  relatively 
complex,  an  approximate  formula  has  been  developed.  This 
approximation  is  called  t^  and  t^  Q°.  Another  approximation 

o 

(z^)  xs  available  such  that  z i  .  However,  z^  is  not  as 

accurate  as  t^ .  Whenever  the  full  computation  is  to  be  omitted, 

it  is  recommended  that  t.  be  used. 

1 

The  observed  data  of  example  1,  part  1,  will  be  used.  Thus, 


a  ~  .80,  a.  ~  .08,  a ..  =  .05,  a  _  =  .01,  a.  *  .005,  a  _ 

0  1  2  3  4  5 


.005 


The  calculations  are  in  three  sections: 


• 

The 

calculation 

of 

fci  > 

• 

The 

calculation 

of 

zi  < 

• 

The 

exact  value 

of 

Q°. 

1.  Calculation  of  t^  (t^_>  ) 

(1)  Calculate  u^.,  the  positive  root  of  equation  50: 


a  a  .  a 

-E  +  -SJi  +  ...  +  n 

u  u^  u 


n-r+1 


1  ”  ao  “  •  *  *  "  ar-l  ’ 


For  r  «  1,  we  obtain 

a  i  a2  a3  a4  a5 

_+  +  +_  + 

u  u  u  u  u 


1  -  ao  ' 
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which  reduces  to 


. 20u5  -  . 08u4  -  . 05u3  -  . Olu2  -  .005u  - 
uj  =  .851  . 


which  reduces  to 


. 12u4  -  .05u3  -  .Olu2  -  . 005u  =  0 
u ^  =  .722  . 

For  r  =  3, 


which  reduces  to 

. 07u3  -  .Olu2  -  . 005u  -  .005  =  0 
u^  =  .531  . 

For  r  =  4, 


which  reduces  to 

. 06u2  -  . 005u  -  .005  =  0 

uj  =  .333  . 

4 

(2)  t^,...,t^  are  given  by  equation  54 


t.  =  Min  ^(u*)1,  (up1-1,...,  ( u|_1 ) 2,-  (u!)j 
We  have 

u!  =  .851,  ul  =  .722,  ul  =  .531,  ul  =  .333  . 
12  3  4 

Hence, 


t2  =  Min  [(up]  =  uj^ 

=  .851 

Min  [{up2,  (up] 

Min  [.724,  .722  ] 

.722 

Min  [(up3,  (up2,  (up] 

Min  [.616,  .521,  .531] 

.521 

Min  [(up4,  (up3,  (up2,  (up] 

Min  [.524,  .376,  .282,  .333]  • 

.282 

t_  is  not  calculated  since  the  exact  value  of 
easily  obtained. 

2.  Calculation  of  (zi  <  Q°) 

The  following  values  must  be  obtained: 
qQ,  the  root  of  equation  26A 


This  has  already  been  obtained  as  uj  .  Thus  qQ  * 
values  of  x^,...,x5  have  been  calculated  in  part  1: 
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Q,-  can  be 


851.  The 


001 


x,  =  .030,  x „  =  .013,  X-  =  .004,  x.  =  .002,  x 
1  2  3  4  5 


+  *2  +  *  '  *  +x^. 


A  ■  xx  -  .030 

A2  =  Xx  +  x2  =  .043 

A^  =  +  X^  +  x^  =  .047 

A^  —  +  x2  +  x^  +  x^  =  .049 

A_  =  x ,  +  x_  +  x,  +  x.  +  xc  *=  .050. 

5  1  2  3  4  5 


From  equation  62  the  lower  bounds  z^  are  calculated: 


zi  "  1  "  1  -  a  -  -  a/;  <  Q?  ‘ 


i-1 


Then 


1  ,  .030  _  _ 

Z1  ~  1  ”  1  -  a  ~  1  ~  .20  ~  * 850 


Z2  =  1  '  1  -  ao  -  ax 
A, 


23  =  1  '  1  -  ac  -  ax  -  a2 

_ N _ 

z4~  1  -  a_  -  a,  -  a„  -  a 

0  12  3 


!  *043  r  A  r\ 

■  1  -  ~rn  *  -642 


=  i  _  • 0^ Z  =  329 

1  .07 


-  1  -  ^  '  -1U3 


zc  is  not  calculated  since  Q°  can  be  obtained  directly, 
5  5 


3.  The  Exact  Value  of  Q7 


We  have  calculated  t^  and  z^  such  that 


-32- 


(i  *  1>  2| • • • r 5) 


Zi  <  Qi  <  *i 
The  exact  value  o£  Q°  is  obtained  as  follows: 

„  u.  I.  f,i-r+l  .  j  *  i-r  o  |  ,  o  .Ti-rl 

"ir"  Min  J(u^)  ,  (u'+1)  ,  ulr  ^c(nlr)|  j  , 

where  u°r  and  ^r(uir)  wiH  be  defined  below. 


°i  - Min  [Hii . 


or  combining  these  equations  with  the  definition  of  t^  we  obtain 


Q°  =  Min  {tj  *  .851 

Q°  =  Min  jt2,  u°21l  ^(u^)]} 

Q°  =  Min  {t3>  n°31l  ^(u^)]2,  u°2[  *2(u°2)]| 

Q°  =  Min  jt4,  u°x[  ^(u^)]3,  u°2[  <j,2(u°2)j2,  u°3[  ♦3(u°3)]} 

If  u°r  >  1,  ^r(u°r)J  >  or  <  4»r(u°r),  then 


u°  <p  (u°  )  1- 
lr  Yr'  ir 


is  neglected  in  the  equations  above. 


!  _  aQ  -  ai  -  a2  -  a3  -  a4 


.005 

.055 


*  .091. 


In  the  equation  of  qV  the  additional  quantities  we  have  to 
compute  are 


♦l(u21) 
,  o  * 
<t»1(u31) 

^2^32* 

^1  ( U  4 1 ) 

VU42} 
***3 ( U  4  3 J 
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The  following  equations  have  exactly  one  positive  root  in  q*,  q*. 
The  root  in  q*  is  u?^;  the  root  in  q*  is  ^(u1?^). 


a*  a*  a 

a?  +  —  +  „  +  . . .  +  - 

q*2  (q$)  <q$> 


7^1  =  (1  '  ao*,qi  ' 


where  q*  satisfies 


(i*  -  l)a*  + 


<i*  -  2)a*  (i*  -  3)a* 

-  +  - 2~ 

q*2  <q$r 


(i*  -  n* ) a** 

+  - =  0, 

(q*,n  _1 


where 


n*  =  n  -  r  +  1 


a* 

o 

=  a  + 
o 

a^  ^  ... 

+  Cl  - 

r-1 

a* 

3 

=  3  j+r 

-1 

(j 

*  i. 

2, . . . , n 

*) 

i* 

=  i  - 

r  +  1  . 

The 

details  of 

the  computation 

are  given  in 

tables 

2  and 

3. 

TABLE  2 

o 

uir 

i 

r 

1  3» 

1 

fo_ 

al 

a* 

3  2 

a3 

a4 

!L 

o 

U21 

2 

1 

5  2 

.80 

.08 

.05 

.01 

.005 

.005 

o 

U31 

3 

1 

5  3 

.80 

.08 

.05 

.01 

.005 

.005 

o 

U32 

3 

2 

4  2 

.88 

.05 

.01 

.005 

.005 

o 

U41 

4 

1 

5  4 

.80 

.08 

.05 

.01 

.005 

.005 

o 

U42 

4 

2 

4  3 

.88 

.05 

.01 

.005 

.005 

o 

U43 

4 

3 

3  2 

.93 

.01 

.005 

.005 

where 

ao  5 

■  .80, 

al  = 

.08,  a2  * 

#05,  3 

3  =  * 

01,  a4 

=  .005, 

a5  = 

.005 

34 


Mutrlctl  Equation  IWiult 


(P*n un°°3)  C  TM1 


Substituting  the  values  from  table  3  in  equation  A  and  neglecting 
several  terms  as  explained  in  table  3,  we  have 


Q°=  .851 

Q°  =  Min  j.  722,  .721} 
Q°  =  Min  {.521,  .517} 

Q°  =  .282 

Q°  =  .091 


.721 

.517 


obtained 

are  shown  in  table 

TABLE  4 

4. 

i 

z . 

1 

0? 

t. 

l 

1 

.851 

.851 

.851 

.851 

2 

.642 

.721 

.722 

.724 

3 

.329 

.517 

.521 

.616 

4 

.183 

.282 

.282 

.524 

5 

.091 

.446 

Thus,  with  the  observed  data,  this  example,  if  all  the 
information  available  about  the  q^'s  is  that 

1  <*2  >  •••  1^5  ' 


all  we  can  say  about  the  is  that 


Q1  1 


.85, 


®2  — 


•  72, 


°3  * 


.52, 


°4  > 


.28,  Q, 


.09 


Note  that 


This  is  always  true 


It  is  interesting  to  compare  with  the  values  of  obtained 
under  the  assumption  that  all  the  q^s  are  equal  and  have  the 
value  qQ.  Under  this  assumption , 

Qj  ■  qQ  ( i  ~  1#2/. ..,5). 

In  table  4,  Q°  =  qQ  and  Q°  is  very  close  to  q^  .  Q°  and  q^ 

differ  by  approximately  .1  and  the  agreement  between  Q°  and 

q*  gets  progressively  worse.  It  will  usually  be  true  that  q* 

and  Q°  are  approximately  equal  for  small  values  of  i,  but  will 
differ  widely  as  i  increases. 


PART  IV 


MINIMUM  AND  MAXIMUM  VALUE  OF  THE  PROBABILITY  THAT  A  PLANE 
WILL  BE  DOWNED  BY  A  GIVEN  NUMBER  OF  HITS  CALCULATED  UNDER 
SOME  FURTHER  RESTRICTIONS  ON  THE 

PROBABILITIES 

In  parts  I,  II,  and  III  we  merely  assumed  that  q^  ,>  q2  —  *  *  *  —  q 

In  many  cases  we  may  have  some  further  a  priori  knowledge 
concerning  the  values  q^,...,qn.  We  shall  consider 

here  the  case  when  it  is  known  a  priori  that  X^qj£  *3j  +  i  <.  *2qj 

(j  =  l,...,n-l),  where  X^  and  X2  ( X^  <  X2  <  1)  are  known 

positive  constants. 

We  shall  also  assume  that 


<  1 


a 


o 


(63) 


Since  aA  +  a2  +  . . .  +  an  <  1  -  aQ,  the  inequality  in  equation 
63  is  certainly  fulfilled  if  X^  is  sufficiently  near  1.  It 
follows  immediately  from  equations  63  and  26  that  q^  <  1. 

CALCULATION  OF  THE  MINIMUM  VALUE  OF  Q.  *  1  -  P.  (i  <  n) 

1  l 

Let  q°,...,q°  be  the  values  of  q^,...,qn  for  which  becomes 
a  minimum.  We  shall  prove  the  following. 

Lemma  1 :  The  relations 

qj+l  =  X2qj  <3  *  (64) 

must  hold. 


Proof:  Suppose  that  the  relation  in  equation  64  does  not  hold 
for  at  least  one  value  j  >  i  and  we  shall  derive  a  contradiction 


J-This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
based  on  Damage  of  Survivors"  was  published  as  SRG  memo  89  and 
AMP  memo  76.4. 
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in  contradiction  to  our  assumption  that  q°  ...  q°  is  a  minimum. 
Hence,  Lemma  1  is  proved. 

Lemma  2 ;  If  j  is  the  smallest  integer  such  that  q£+i  =  *2qk  for 
all  k  >  j,  then  q°  ■  A^q®^  for  r  *  2,3,...,j-l. 


Proof:  Assume  that  Lemma  2  does  not  hold  and  we  shall  derive  a 
contradiction.  Let  u  be  the  smallest  integer  greater  than  one 

such  that  q°  >  A^q®^  •  It  follows  from  the  definition  of  the 

integer  u  that  if  u  >  2,  then  q°_^  =  From  assumption  63 

it  follows  that  q°  <  1.  Hence,  if  we  replace  q°_^  by 

q^_^“  (1  +  e)q°_i  (e  >0),  then  for  sufficiently  small  e  the 
inequalities  Ajqr  £  qr+1  <.  A2qr  (r  =  l,...,n-l)  will  not  be 
disturbed.  Let  v  be  the  smallest  integer  greater  than  or  equal 
to  u  such  that  q°+i  <  A2q°  •  Since  by  assumption  j  is  the 


smallest  integer  such  that  q£+i  *  A2q£  for  k  1.  3»  WG  must 


have  q j  <  A2qj_j_.  Hence,  v  <  j-1.  It  is  clear  that  replacing 
o  q° 

qy  by  q^  ■  ^  ^  —  we  shall  not  disturb  the  inequalities 
Ajqr  £  qr+i  £  A2q^,  ^  *  !»•••» n—1 ) .  Hence ,  if 
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o 

% 

.  +  E 


qu-l  -  (1  +  E)  qu.1#  q'  =  ,  and  q* 


for  r  /  u,  /  v,  then  X^  <  q£+1  <  *2q£  (k  =  l,...,n-l)  is  ful¬ 
filled.  Furthermore,  we  have 


qi 


•  o 

q[  -  qA 


o  , 
q^  and 


n 

E 

j=i 


qi 


qj 


<  1  -  a 


Hence,  there  exists  a  positive  A  <  1  such  that 


n 


E 

j=i 


i 


-  a 


o 


and  q!j  ■  Aq^  (j  =  l,...,n).  But  then 

qj  ...  q£  <  q{  ...  q-  -  q°  •••  q° 

in  contradiction  to  the  assumption  that  q°  ...  q°  is  a  minimum. 

Hence,  Lemma  2  is  proved. 

Let  E.  (r  =  l,...,i-l)  be  the  minimum  value  of  Q.  under  the 
ir  i 

restriction  that  q^+i  =  ^2qj  for  ^  “  r+l,...,n-l  and  q^+i  = 
for  3  =  l,...-,r-l.  From  Lemma  1  and  2  it  follows  that  the  mini¬ 
mum  of  Q.  is  equal  to  the  smallest  of  the  i  -  1  values  E. ,,..., E.  .  , 

X  11  1 | 1~1 

The  computation  of  the  exact  value  of  Eir  can  be  carried  out 

in  a  way  similar  to  the  computation  of  Mir  described  in  part 

II.  Since  these  computations  are  involved  if  n  is  large,  we  shall 
discuss  here  an  approximation  method. 

Let  Ejr  (r  *  l,...,i-l)  be  the  value  of  if  q^  +  1  *  *2qj  for 
j  *  r+l,...,n-l  and  q^+1  *  for  j  *  l,...,r.  Furthermore, 

let  e£q  be  the  value  of  if  q^  +  1  ■  X2q^  (j  ■  l,...,n-l).  Then, 

if  n  is  large,  the  minimum  of  Ef  _  .  and  E*  will  be  nearly  equal 
to  E^.  Hence,  we  obtain  an  approximation  to  the  minimum  of 
by  taking  the  minimum  of  the  i  numbers  E*Q,  EJ^,...,EJ 
The  quantity  E^r  can  be  computed  as  follows.  Let  gr  be  the 
positive  root  in  q  of  the  equation 


r+1 

V* 

a.  n-r-1 

1  i  V 

ar+l+j 

L. 

j=l 

ft 

»1  o3 

r(r+l)  j  j(j+l) 

»x  2  3  X2  2  q"1*2 

(r  =  0, 1, . . . , i- 

Then 

E*  = 
lr 

£if±ILcti-r-l) 
X1  2 

(i-r)(i-r-l) 

*  2  i 

\2  gr  • 

(67) 


MINIMUM  OF  Q 

n 


Let  q?,...,q°  be  values  of  q,,...,q..  for  which  Q„  becomes  a 
minimum.  We  shall  prove  that  q°+i  =  *j_q°  (j  =  l#*».»n-l). 

Assume  that  there  exists  a  value  j  <  n  such  that  q°+^  >  ljq° 

and  we  shall  derive  a  contradiction.  Let  u  be  the  smallest 
integer  such  that  q°+1  >  *jq°  and  let  v  be  the  largest  integer 

q° 

such  that  q°+1  >  1^°  .  Let  q^  =  (1  +  e)q°  (e  >0),  q^+1=  T^Tc' 

and  q!  =  q°  for  j  ^  u,  /  v+1.  Then  for  sufficiently 
small  e  we  shall  have  A^q^.  <  *1^+1  —  ^2^r  =  • • • rn_l) • 

Furthermore,  we  have 


*1 


q°  and 
n 


<  1 


Hence,  there  exists  a  positive  A  <  1  such  that  q? 
(j  =  l,...,n)  and 


1 


-  a 


o* 


qj 
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But  then  <  q°  . . .  q°  in  contradiction  to  the  assumption 

that  q°  ...  q°  is  a  minimum.  Hence,  our  statement  is  proved. 

If  q  is  the  root  of  the  equation 


n  a . 


n(n-l) 

then  the  minimum  of  Qn  is  equal  to  X^  2  qn 
MAXIMUM  OF  Q.  (i  <  n) 

Let  q*,...,q*  *3e  values  of  q^».../qn  for  which  becomes  a 
maximum.  We  shall  prove  the  following! 

Lemma  3:  The  relations 


qj+l  =  Xlqj  (j  * 

must  hold. 

Proof:  Assume  that  there  exists  an  integer  j  _>  ■*■  suc^  that 
+  1  >  X^q*  and  we  shall  derive  a  contradiction.  Let  qjj.  «*  q 

for  r  =  l,...,i  and  let  q^  +  1  =  X^q.!  (j  =  i,...,n-l).  Then 


q{  •••  q|  -  qj  • 


j?x  qi  •••  qj 


>  1  -  a. 


Hence,  there  exists  a  value  A  >  1  such  that 


1 


9 


where  q!j  ■  Aq^  (j  =  l,...,n).  But  then  qj  ...  q?  >  q|  ...  q£  in 
contradiction  to  the  assumption  that  q*  ...  q*  is  a  maximum. 
Hence,  Lemma  3  is  proved. 
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fi » 


Lemma  4 : 


If  for  some  j  <  i  we  have  q^  +  1  >  X^q * ,  then 
qk+l  =  X2qk  for  k  = 

Proof:  Assume  that  q* , ,  >  X.q*  for  some  j  <  i  and  that  there 

^3+1  I]  j 

exists  an  integer  k  £  j-1  such  that  q£  +  1  <  X2q*  *  We  shail 

derive  a  contradiction  from  this  assumption.  Let  u  be  the 
smallest  integer  such  that  q*  +  ^  <  X2q*  •  Furthermore,  let  v  be 

the  smallest  integer  greater  than  or  equal  to  u  +  1  such  that 

q* 

qv+l  >  Xlqv  *  Ifc  is  clear  that  v  <.  3*  Let  ^  *  TTe  (e  >  °>  • 
q^  =  (1  +  e)  q*  ,  and  qj_  =  q*  for  r  /  u,  /  v.  Then  for  suffi¬ 
ciently  small  e  we  have 

Xlqj  1  qj+l  i  X2qj  {j  =  1'***'n“1)  * 


Furthermore,  we  have 


qi 


qi 


=  a* 


qi 


and 


n 

Z.  qj- 

D  =1  M1 


’j 


>  1  -  a. 


Hence,  there  exists  a  value  A  >  1  such  that 


j=i q: 


a . 


qj 


=  1  -  a. 


where  q?  =  Aq!  (j  =  i,...,n). 


But  then  q' 


contradiction  to  the  assumption  that  q| 


Ji 

is 


>  • 
a  maximum. 


qi 


in 


Let  (r  »  l,...,i-l)  be  the  maximum  of  under  the  restric¬ 

tion  that  q.  ,  =  X,q.  for  j  =  r+l,...,n-l  and  q.  ,  =  X„,q  .  for 
D+l  1  ]  j+1  2  ] 

j  =  l,...,r-l.  From  Lemma  3  and  4  it  follows  that  the  maximum  of 
Qi  is  equal  to  the  maximum  of  the  i  -  1  values  D^1,...,Di 

The  computation  of  the  exact  value  of  can  be  carried  out  in  a 

way  similar  to  the  computation  of  M^r  in  part  II.  Since  these 

computations  are  involved  if  n  is  large,  we  shall  discuss  here 
only  an  approximation  method. 
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Let  D£r  (r  *  l,...,i-l)  be  the  value  of  if  c3j  +  ^  ■  X^q^  ^or 

j  *  r+l,...,n-l  and  q^+1  =  X2q^  for  j  *  l,...,r.  Furthermore, 

let  Djo  be  the  value  of  if  q..  +  1  =  X^q..  (j  *  l,...,n-l).  Then, 

if  X.  is  not  much  below  one,  the  maximum  of  D*  and  D*  . 
x  ir  if r* x 

(r  =  l,...,i-l)  will  be  nearly  equal  to  D^r.  Hence,  we  obtain  an 

approximation  to  the  maximum  value  of  Q.  by  taking  the  largest  of 

the  i  values  D?  , ...,D*  .  .. 

1 0  1 1 1“1 

The  value  of  D*r  can  be  determined  as  follows.  Let  gr  be  the 
root  in  q  of  the  equation 


r+1 


a . 


n-r-1 


&  JUjll.  + 

1  '  r»J 


r+l+j 


0  ■+jr 


=  1  -  a. 


2  r+l+j 


Then 


D* 

ir 


^f^  +  U-r-l)  Ii.r£-1)  (i-D 

4  xi  K 


MAXIMUM  OF  Qn 

We  shall  prove  that  the  maximum  of  Qn  is  reached  when  qj+1=  X2qj 

(j  =  l,...,n-l).  Denote  by  qj  ...  q*  the  values  of  q^  ...  qn 

for  which  Qn  becomes  a  maximum.  We  shall  assume  that  there 

exists  a  value  j  <  n  such  that  q*+1  <  X2q*  and  we  shall  derive 

a  contradiction  from  this  assumption.  Let  u  be  the  smallest  and 
v  be  the  largest  integer  such  that  q*  +  1  <  X2q*  and  q*+1  <  X2q*  . 

q* 

Let  q^  =  l  +Ue  {£  >0),  q^+1  -  { 1  +  e)  q£+1,  and  q£  =  q*  for 
r  /  u,  /  v+1.  Then  for  sufficiently  small  e  we  shall  have 
Xlqr  ^  qr+l  i  X2qr  (r  *  ^•••»n-l). 
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Furthermore,  we  have 


qi  qA  =  q!  • 


.  q*  and  £ 

j=l 


n  \  ^  •  q 


7  >  1  “  a. 


Hence,  there  exists  a  value  A  >  1  such  that  qV  =  Aq^ 
( j  =  1, . . . , n )  and 


n 

Z>.  77" 


a  j 


j— 1  Mj 


tt  =  1  -  a  . 
q"  o 


But  then  q£  ...  q^  >  q*  ...  q*  in  contradiction  to  the  assumption 
that  q*  ...  q*  is  a  maximum.  Hence,  our  statement  is  proved. 

The  maximum  of  Qn  is  equal  to 


n(n-l) 


where  q  is  the  root  of  the  equation 


NUMERICAL  EXAMPLE 

The  same  notation  will  be  used  as  in  the  previous  numerical 
examples.  The  assumption  of  no  sampling  error,  which  is  common 
to  all  the  previous  examples,  is  retained.  In  part  I  it  was 
assumed  that  the  q^ ,  the  probability  of  a  plane  surviving  the 

i-th  hit,  knowing  that  the  first  i  -  1  hits  did  not  down  the 
plane,  were  equal  for  all  i  (q^  =  q2  =  . . .  =  qn  =  qQ  (say)). 

Under  this  assumption,  the  exact  value  of  the  probability  of  a 
plane  surviving  i  hits  is  given  by 


In  part  III  it  was  assumed  that  q.  >  q„  >  ...  >  q  .  Since  no 
lower  limit  is  assumed  in  the  decrease  from  q.  to  q . , , ,  only  a 

l  ^i+l  1 
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lower  bound  to  the  could  be  obtained.  The  assumption  here  is 
that  the  decrease  from  to  q^+^  lies  between  definite  limits. 
Therefore,  both  an  upper  and  lower  bound  for  the  can  be 
obtained . 

We  assume  that 


Vi  i  qi+l  <  Vi  ’ 

where  X^  <  X2  <  1  and  such  that  the  expression 


n 

£ 


JL 


j(j-D 


<  1  -  a 


(A) 


is  satisfied. 

The  exact  solution  is  tedious  but  close  approximations  to  the 
upper  and  lower  bounds  to  the  0^  for  i  <  n  can  be  obtained  by 

the  following  procedure.  The  set  of  hypothetical  data  used  is 


a_  = 


a,  ■ 


a„  = 


X.  = 


.780 

.070 

.040 

.80 


a4  = 
a5  = 
X2 


.010 

.005 

.005 

.90 


Condition  A  is  satisfied,  since  by  substitution 


.07 


+ 


.04  .01  +  .005  +  .005 

.8  (.8)3  (.8)6  ( . 8 ) 10 


which  is  less  than 


.20529  , 


1  -  a„  =  .22  . 
o 


THE  LOWER  LIMIT  OF  Qi 

The  first  step  is  to  solve  equation  66.  This  involves  the 
solution  of  the  following  four  equations  for  positive  roots  g^, 


4  7- 


;  ■  n  a  itiSva&iMiMXtk  • » 


91'  g2'  9 3 


Ik  +  _!l_  +  43  +  ®4  +  — f  5  _  =  !  _  a 

q  Xiq2  Xjq4  xjx2q5  C 


.07  .  .04  .  .01  .005  .005  _  _  _  ,, 

~  — ■  +  “■  .  +  '  1  v  +  1  +  c  “  •  ^  ^ 

q  .Bq^  . 512q J  .262144q*  (. 134218) (. 9 )q3 

. 22q5  -  .07q4  -  .05q3  -  .019531q2  -  .019073q  -  .041392  =  0 


g^  =  .964  . 

Next,  calculate  the  i  numbers  defined  by 


E*  =  xa(i,r)  xb(i'r)  g1 

Eir  A1  A2  yr 


(r  «  0,1, ... , i-1) , 


where 


a ( i ,r )  -  r(y  1-)-  +  r(i  -  r 
b  ( i  ,r )  g  (i  -  r)(i  -  r  -.1) 


.844 

.904 

.941 

.964 


The  minimum  of  the  E?r  (r  »  0,...,i-l)  will  be  the  lower  limit  of 
The  computations  are  given  in  table  5. 


The  lower  limit  of  Q5  can  be  obtained  directly.  The  lower  limit 
of 


°5  "  *i0c35  * 


where  q  is  the  positive  root  of 


al  a2  a3  a4  a5 
“  +  ,  2  +  7T“5  +  76~4  +  >10  5 

q  xxq  ^q  xxq  xA  q 


1  -  a. 


.07  .04  +  .01  +  .005  +  .005 

q  . 8q2  . 512q3  .262144q4  .107374q5 


.22 


q  =  .974  • 

The  lower  limit  of 

Q5  *  ( . 8) 10( .974) 5  *  .094  . 

THE  UPPER  LIMIT  OF 

The  computations  for  the  upper  limit  of  are  entirely  analogous 

to  the  computations  of  the  lower  limit.  First ,  we  solve  the 
equations  of  part  IV,  which  for  this  example  are  the  following: 


al  a2  .  a3  a4  .  a5 

_  +  ,  2  +  ,3  3  +  ,6  4  +",i075  85  1  “  ao 

q  M  x,q  x,q  x,  q 


.07  .  .04  .  .01 

"  "  T  a  +  -i  T 


.005 


.005 


q  .8q2  .512q3  .262144q4  ,107374q5 


.22 


. 22q5  -  . 07q4  -  .05q3  -  ,019531q2  -  ,019073q  -  .046566  »  0 


.974 
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.07  .04  .01  .005  _ .005 _  =  22 

q  . 9q2  .729q3  (. 59049) ( -8)q4  ( . 512) ( .478297)q5 

. 22q5  -  . 07q4  -  .044444q3  -  .013717q2  -  .010584q  -  .020417  =  0 
g*  =  .869 


a 


o 


.07  .04  .01  .005  .005 

-  +  - r  +  - r  +  - t  +  - c  =  -22 

q  . 9q^  . 7  2  9q  J  .531441q  (. 387420 )(. 8 )q 

. 22q5  -  .07q4  -  .044444q3  -  .0137l7q2  -  .009408q  -  .016132  =  0 
g*  »  .851 
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Next,  calculate  the  i  numbers  defined  by 


a(i,r)  b(i,r)  *i 
2  X1  gr 


(r  -  0,l,...,i-l), 


where 


a(i,r)  »  r-(‘r2+  ^  +  r(i  -  r  -  1) 


h/i  y\  -  (i  ~  r><*  ~  r  ~  1) 

b(i,r)  —  —  r  _  2 


g*  -  .974 
g*  *  .905 
g*  *  .869 
g*  »  .851 

The  maximum  of  the  D*r  (r  =  0,...,i-l)  will  be  the  upper  limit  of 

Qi»  The  computations  are  given  in  table  6. 

The  upper  limit  of  Q5  can  be  obtained  directly.  The  limit  of 


The  upper  limit  of 

Qs  -  ( . 9 ) 1 0 C .844) 5  -  .149 

Summarizing  the  results,  the  upper  and  lower  limits  of  the 
probability  of  a  plane  surviving  i  hits  are  given  by 

.844  <  Q  <  .974 
.641  <  Q2  <  .759 
.426  <  Q3  <  .480 

.226  <  Q.  <  .279 

4 

.094  <  Qg  <  .149 
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PART  V 


SUBDIVISION  OF  THE  PLANE  INTO  SEVERAL 
EQUI- VULNERABILITY  AREAS1 

In  parts  I  through  IV  we  have  considered  the  probability  that  a 
plane  will  be  downed  by  a  hit  without  any  reference  to  the  part 
of  the  plane  that  receives  the  hit.  Undoubtedly,  the  probability 
of  downing  a  plane  by  a  hit  will  depend  considerably  upon  the 
part  that  receives  the  hit.  The  purpose  of  this  memorandum  is  to 
extend  the  previous  results  to  the  more  general  case  where  the 
probability  of  downing  a  plane  by  a  hit  depends  on  the  part  of 
the  plane  sustaining  the  hit.  To  carry  out  this  generalization 
of  the  theory,  we  shall  subdivide  the  plane  into  k  equi- 
vulnerability  areas  A. , . . . , Ak .  For  any  set  of  non-negative 

integers  i^,...,i^  let  P( i^ , . . . , i^)  be  the  probability  that  a 
plane  will  be  downed  if  the  area  A^  receives  i^  hits,  the  area 
A2  receives  i2  hits,...,  and  the  area  A^  receives  i^  hits.  Let 
Q(i1,...,ik)  =  1  -  P(i1,...,ik).  Then  Q(i1,...,ik)  is  the  prob¬ 
ability  that  the  plane  will  not  be  downed  if  the  areas  A^,...,Ak 
receive  i^,...,i  hits,  respectively.  We  shall  assume  that 
Q(i1,...,ik)  is  a  symmetric  function  of  the  arguments  i^,...,ik. 


To  estimate  the  value  of  Q(i^,...,ik)  from  the  damage  to 

returning  planes,  we  need  to  know  the  probability  distribution  of 
hits  over  the  k  areas  A^,...,Ak  knowing  merely  the  total  number 

of  hits  received.  In  other  words,  for  any  positive  integer  i  we 
need  to  know  the  conditional  probability  ( i^ , . . . , ik )  that  the 

areas  A1#...,Ak  will  receive  hits,  respectively, 

knowing  that  the  total  number  of  hits  is  i.  Of  course, 

( ii , . . . , ik)  is  defined  only  for  values  i^,...,ik  for  which 

i^  +  ...  +  ik  =  i  .  To  avoid  confusion,  it  should  be  emphasized 

that  the  probability  Y^ ( , . . . r i^)  is  determined  under  the 


1This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  96  and 
AMP  memo  76.5. 


assumption  that  dummy  bullets  are  used.  It  can  easily  be  shown 
that  it  is  impossible  to  estimate  both  \  ( i ^ > • . . , i^ )  and 

Q( i^ , . . . , i^)  from  the  damage  to  returning  planes  only.  To  see 

this,  assume  that  k  is  equal  to  2  and  all  hits  on  the  returning 
planes  were  located  in  the  area  A^ .  This  fact  could  be  explained 

in  two  different  ways.  One  explanation  could  be  that 

Yf ( i i # i 2 )  =  0  for  i2  >  0.  The  other  possible  explanation  would  be  that 

Q(il,i2)  *  0  for  i2  >  0.  Hence,  it  is  impossible  to  estimate 

both  Yi(i1,i2)  and  Q(ii,i2)*  Fortunately,  Y^ij^. . . ,  i^)  can  be 

assumed  to  be  known  a  priori  (on  the  basis  of  the  dispersion  of 
the  guns),  or  can  be  established  experimentally  by  firing  with 
dummy  bullets  and  recording  the  hits  scored.  Thus,  in  what 
follows  we  shall  assume  that  Y^ ( i^ , . . . , i^ )  is  known  for  any  set 

of  integers  i.,,..,ik. 


Clearly,  the  probability  that  i  hits  will  not  down  the  plane  is 
given  by 

3  S  •••  51  ( i^,  •  •  • ,  ijj  )Q(  i^ » •  •  • » ijj )  #  (69) 

lk  11 

where  the  summation  is  to  be  taken  over  all  non-negative  integers 
i^,...,^  for  which  i^  +  . . .  +  i^  =  i  . 


Let  6^  (ijy...,ik)  be  the  conditional  probability  that  the  areas 

A^,...,Ak  received  i^,...,^  hits,  respectively,  knowing  that  the 

plane  received  i  hits  and  that  the  plane  was  not  downed.  Then  we 
have 


fii^*l"  *  * ' *k ^ 


Y^(i^,.*.,ij^)Q(i^,...,ij^) 

*1 


(70) 


Of  course,  6^ ( i ^ , . . . , i^)  is  defined  only  for  non-negative 
integers  for  which  iA  +  . . .  +  i  ■  i  . 
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The  probability  6 . ( i  ,  ...,i  )  can  be  determined  from  the  distri- 

bution  of  hits  on  returning  planes.  In  fact,  let  a( i^ . i^) 

be  the  proportion  of  planes  (out  of  the  total  number  of  planes 
participating  in  combat)  that  returned  with  i^  hits  on  area  , 

i 2  hits  on  area  • • • ,  and  i^  hits  on  area  A^.  Then  we 

obviously  have 


6i(i1, . . ,,xk)  - 


t  •  •  •  t  i]^ ) 


From  equations  70  and  71,  we  obtain 

.  %  _  Qia(ii . ik) 

QU1 . 1k)  "  a.Y.TT".-:.,'^) 


(i  =  i,  +  . . .  +  ik) 


Since  can  be  estimated  by  methods  described  in  parts  I  through 
IV,  estimates  of  Q(i^,...,ik)  can  be  obtained  from  equation  72. 

According  to  equation  29,  the  probabilities  Q^,...,Q  satisfy  the 
equation 


n  a . 

T  ~2  =  i  _  a 

A.  Q  .  c 

3=1  3 


We  have  assumed  that  q^^  _> 

stating  that 


.  qn .  This  is  equivalent  to 


Qi+1 


Q.  -  Q  ■ 
i  3 


for  j  <  i 


A  similar  assumption  can  be  made  with  respect  to  the  prob¬ 
abilities  0(1^,  . .  .,±k)  .  In  fact,  the  conditional  probability 

that  an  additional  hit  on  the  area  Af  will  not  down  the  plane 

knowing  that  the  areas  A^ , . . . , A^  have  already  sustained 

i^,...,ik  hits,  respectively,  is  given  by 


-58- 


(75) 


Q(i1,...#ir_1/  ir+l,  ir+1, 


V 


Obviously#  we  can  assume  that  if 


jl-  1l'  j2  -  i2'***'jk  -  i) 


then 


0(  11 # . . • # iv^ # ir+l, ir+1 , . . . , ik )  Q( jr,1#jr+l#jr+I#...#jk) 

Q(i1,..."#ir_1,ir,ir+1,.../ik)  -  "  Q( j1,...,3r_1»3r#3r+1#***#3k) 


for  r  =  1, 2, . . . ,k. 


(76) 


Hence#  the  possible  values  of  Q1#...,Qn  are  restricted  to  those 

for  which  equation  73  is  fulfilled  and  for  which  the  quantities 
Q(i^#...#ik)  computed  from  equation  72  are  less  than  or  equal  to 

one  and  satisfy  the  inequalities  of  equation  76.  It  should  be 
remarked  that  the  inequalities  of  equation  76  do  not  follow  from 
the  inequalities  of  equation  74.  From  equation  72  and  the 
inequality  Q(i1#...,ik)  <  1,  it  follows  that 


Q,  < 


ai  V1!'*  *  *  #1k} 


i  a( i^, . . . , ik) 


(77) 


If  the  right-hand  side  expression  in  equation  77  happens  to  be 
less  than  one#  then  equation  77  imposes  a  restriction  on  . 

Since 


E 

4k 


E 

ii 


a(i,  #. . 


Yi(  il#  •  *  •  #  ik) 


1 


(the  summation  is  taken  over  all  values  i^ , . . . , ik  for  which 
ij+  ...  +  ik  ■  i)#  we  must  have  either 


i 
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1 


ai  w* 


"V 


for  all  values  i^,...,ik  for  which  +  ...  +  ik  =  i,  or 


aiYi(il 


a  ( l 


.  #  iR) 


ik) 


1 


at  least  for  one  set  of  values  , . . . , ik  satisfying  the  condition 
i  +  ...  +  i.  =  i  .  Hence,  equation  77  gives  an  upper  bound  for 
Qi  whenever  there  exists  a  set  of  integers  i^,...,ik  such  that 
i^  +  . . .  +  ik  =  i  and 


a(i. 


9  • 


"1k) 


^Y,  (i 


It  is  of  interest  to  investigate  the  case  of  independence,  i.e., 
the  case  when  the  probability  that  an  additional  hit  will  not 
down  the  plane  does  not  depend  on  the  number  and  distribution  of 
hits  already  received.  Denote  by  q(i)  the  probability  that  a 
single  hit  on  the  area  A^  will  not  down  the  plane.  Then  under 

the  assumption  of  independence  we  have 

Q(i1,...,ik)  =  [q( 1 ) ] 11  [q{ 2 ) ] 1 2  ...  [q(k)]  *  .  (78) 


Hence,  the  only  unknown  probabilities  are  q( 1 ) , . . . ,q(k) . 

Let  Y(i)  be  the  conditional  probability  that  the  area  A^  is  hit 
knowing  that  the  plane  received  exactly  one  hit.  Obviously 


Y-i  ( i 


'V  = 


il 


1ll 


V 


f Y(l) J 


[ Y(k) ] 


(79) 
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V 


Similarly,  let  5(i)  be  the  conditional  probability  that  the  area 
is  hit  knowing  that  the  plane  received  exactly  one  hit  and 

this  hit  did  not  down  the  plane.  Because  of  the  assumption  of 
independence,  we  have 


*i<  il- 


,1k)  * 


xk! 


[6(1)] 


I6(k)) 


(CO) 


Furthermore,  we  have 


_  Y(i)q(i) 
k 

£  Y(i)q(i) 
i  =  l 


(Cl) 


Since  the  probability  q  that  a  single  hit  does  not  down  the  plane 
k 

is  equal  to  £  Y(i)q(i),  we  obtain  from  equation  81 

i  =  l 


q(  i) 


(  02 ) 


Because  of  the  assumption  of  independence,  we  see  that  6 ( i )  is 
equal  to  the  ratio  of  the  total  number  of  hits  in  the  area  A  of 

the  returning  planes  to  the  total  number  of  hits  receive  by  the 
returning  planes.  That  is 


7!  •••  ^  j  ^ a  (  j  ^  •  •  •  •  •  j  ^ ) 

Jk  J1 

"  ■■  ■  ■  ■  * 

T!  •  •  •  5^  ( j i  •••  + 

Jk  jl 


Since  Y(i)  is  assumed  to  be  known  and  since  $(i)  can  be  computed 
from  equation  83,  we  see  from  equation  82  that  q(i)  can  be 
determined  as  soon  as  the  value  of  q  is  known.  The  value  of  q 
can  be  obtained  by  solving  the  equation 
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NUMERICAL  EXAMPLE 


In  the  examples  for  parts  I,  III,  and  IV  we  have  estimated  the 
probability  that  a  plane  will  be  downed  without  reference  to  the 
part  of  the  plane  that  receives  the  hit.  However,  the  vulner¬ 
ability  of  a  particular  part  (say  the  motors)  may  be  of  interest 
and  this  example  illustrates  the  methods  of  estimating  part 
vulnerabilities  under  the  following  assumptions: 

•  The  number  of  planes  participating  in  combat  is  large  so 
that  sampling  errors  can  be  neglected. 

•  The  probability  that  a  hit  will  down  the  plane  does  not 
depend  on  the  number  of  previous  non-destructive  hits.  That 
is,  qx  =  q2  =  ...  =  qn  =  qQ  . 

•  Given  that  a  shot  has  hit  the  plane,  the  probability  that 
it  hit  a  particular  part  is  assumed  to  be  known.  In  this 
example  it  is  put  equal  to  the  ratio  of  the  area  of  this 
part  to  the  total  surface  area  of  the  plane. ^ 

•  The  division  of  the  plane  into  several  parts  is  repre¬ 
sentative  of  all  the  planes  of  the  mission.  If  the  types  of 
planes  are  radically  different  so  that  no  representative 
division  is  possible,  we  may  consider  the  different  classes 
of  planes  separately. 

Consider  the  following  example.  Of  400  planes  on  a  bombing 
mission,  359  return.  Of  these,  240  were  not  hit,  68  had  one  hit, 
29  had  two  hits,  12  had  three  hits,  and  10  had  four  hits. 
Following  the  example  in  part  I  we  have 

N  =  400, 

whence 


240 

a  =  . 600 

o 

68 

ax  =  .170 

29 

a2  =  .072 

12 

a^  =  .030 

10 

a.  =  .025 
4 

*By  area  is  meant  here  the  component  of  the  area  perpendic¬ 
ular  to  the  direction  of  the  enemy  attack.  If  this  direction 
varies  during  the  combat,  some  proper  average  direction  may  be 
taken . 


As  before,  the 
plane  is  given 


probability  that  a  single  hit  will  not  down  the 
by  the  root  of 


i 


which  reduces  to 


.4q*  -  .  170q^  -  .  072q^  -  . 030qQ  -  .025  =  0 

and 


q  =  . 850  • 
o 

Suppose  that  we  are  interested  in  estimating  the  vulnerability  of 
the  engines,  the  fuselage,  and  the  fuel  system.  Assume  that  the 


following  data 
mission: 

is  representative  of 

all 

the 

planes 

of  the 

Part  number 

Description 

Area  of 

part 

Ratio  of 
area  of  part 
to  total 
area  (Y(i)) 

1 

2  engines 

35 

sq . 

ft. 

35 

130  “ 

.  269 

2 

Fuselage 

45 

sq . 

ft. 

45 

130  “ 

.  346 

3 

Fuel  system 

20 

sq . 

ft. 

20 

130  “ 

.  154 

4 

All  other  parts 

30 

sq . 

ft. 

30 

no  ~ 

.2  31 

Total  area 

130 

sq . 

ft. 

The  ratio  of  the  area  of  the  i-th  part  to  the  total  area  is 
designated  Y(i).  Given  that  the  plane  is  hit,  by  the  third 
assumption,  Y(i)  is  the  probability  that  this  hit  occurred  on 
part  i .  Thus 
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Y(l)  =  .269 
Y(  2 )  =  .346 
Y(  3)  =  .154 
Y (  4 )  =  .231 

The  only  additional  information  we  require  is  the  number  of  hits 
on  each  part.  Let  the  observed  number  of  hits  be  202.  In 
general,  the  total  number  of  hits  (on  returning  planes)  must  be 
equal  to 


A^  +  2A2  +  3A  ^  + 
and  in  this  example 


+  nA 

n 


Al  +  2A2  +  3A3 


+  4A4  =  68  +  2(29)  +  3(12)  +  4(10) 


202 


The  hits  on  the  returning  planes  were  distributed  as  follows: 

Ratio  of  number  of  hits 
observed  on  part  to 

Number  of  hits  total  number  of  observed 

Part  number  observed  on  part  hits  ( 6 ( i )  ) _ 


1 

39 

.193 

2 

78 

.  386 

3 

31 

.154 

4 

54 

.267 

Total  number  of  hits 

202 

The  ratio  of  the  number  of  hits  on  part  i  to  the  total  number  of 
hits  on  surviving  planes  is  designated  5(i).  Then  q(i),  the 
probability  that  a  hit  on  the  i-th  part  does  not  down  the  plane, 
is  given  by 


’,i>  “IflT 
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whence 


_  6(1) 
q(1)  -  wry 

qo  = 

.193 

.269 

(.850)  = 

.61 

_  4(2 ) 
q(2)  -  ITT) 

qo  " 

.386 
.  346 

(.850)  = 

.95 

_  6(3) 

q(  3)  ttit 

qo 

.154 

.154 

(.850)  = 

.85 

q(4)  “  TUT 

qo  = 

.267 
.2  31 

(.850)  = 

.  98 

The  results  may  be  summarized  as  follows: 


Part 


Probability  of 
surviving  a  single 
hit  (q( i ) ) _ 


Probability  of  being 
downed  by  a  single 
hit  (1  -  g(i)  ) _ 


Entire  plane 

.85 

.15 

Engines 

.61 

.  39 

Fuselage 

.95 

.05 

Fuel  system 

.85 

.15 

Other  parts 

.98 

.02 

Thus,  for  the  observed  data  of  this  hypothetical  example,  the 
engine  area  is  the  most  vulnerable  in  the  sense  that  a  hit  there 
is  most  likely  to  down  the  plane.  The  fuselage  has  a  relatively 
low  vulnerability. 
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PART  VI 


SAMPLING  ERRORS -l 

In  parts  I  through  V  we  have  assumed  that  the  total  number  of 
planes  participating  in  combat  is  so  large  that  sampling  errors 
can  be  neglected  altogether.  However,  in  practice  N  is  not 
excessively  large  and  therefore  it  is  desirable  to  take  sampling 
errors  into  account.  We  shall  deal  here  with  the  case  when 
q ^  =  q ^ •  =  qn  =  q  (say)  and  we  shall  derive  confidence  limits  tor 

the  unknown  probability  q. 

If  there  were  no  sampling  errors,  then  we  would  have 

(155) 

xi  =  a0  “  al  •**  ~  ai-l  X1  ”  x2  ***  ”  xi-l^ 


(i  =  2,3,...  )  , 


where  p  =  1  -  q.  However,  because  of  sampling  errors  we  shall 
have  the  equation 


"  ^i(1  - 


a 


-  a 


i-1 


-  x 


-  xi-l>' 


(06) 


where  p^  is  distributed  like  the  success  ratio  in  a  sequence  of 

tk  =  U ( 1  -  aQ  -  a^  -  ...  -  a^_^  -  x  ^  -  ...  -  x^_^)  independent 

trials,  the  probability  of  success  in  a  single  trial  being  equal 
to  p. 

Let  q^  =  1  -  p^  .  Then,  according  to  equation  26  we  have 


z  - 


=  1  -  a  . 


( « >  7  ) 


j=l  q1  ...  q j 


J-This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  103  and 
AMP  memo  76.6. 
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provided  that  x^  =  0  for  i  >  n.  In  part  I  we  have  shown  that 

x^  =  0  for  i  >  n  if  there  are  no  sampling  errors.  This  is  not 

necessarily  true  if  sampling  errors  are  taken  into  account.  However, 
in  the  case  of  independence,  i.e.,  when  q^  =  q  ( i  =  1,2,...),  x^ 

00 

is  very  small  for  i  >  n  so  that  £  x,  can  be  neglected. 

i=n+l 

In  fact,  if  the  number  of  planes  that  received  more  than  n  hits 
were  not  negligibly  small,  it  follows  from  the  assumption  of 
independence  that  the  probability  is  very  high  that  at  least  some 
of  these  planes  would  return.  Since  no  plane  returned  with  more 

oo 

than  n  hits,  for  practical  purposes  we  may  assume  that  ^  xi=0. 

In  what  follows  we  shall  make  this  assumption.  i=n+l 

Each  of  the  quantities  q^ , . . . ,qn  can  be  considered  as  a  sample 

estimate  of  the  unknown  probability  q.  However,  the  quantities 
q^,...,qn  are  unknown.  It  is  merely  known  that  they  satisfy 

the  relation  in  equation  87.  Confidence  limits  for  q  may  be 
derived  on  the  basis  of  equation  87.  However,  we  shall  use 
another  more  direct  approach. 

To  derive  confidence  limits  for  the  unknown  probability  q  we  shall 
consider  the  hypothetical  proportion  b^  of  planes  that  would  have 

been  hit  exactly  i  times  if  dummy  bullets  would  have  been  used. 

We  shall  treat  the  quantities  b^,...,b^  as  fixed  (but  unknown) 

constants,.  This  assumption  docs  not  involve  any  loss  of 
generality,  since  the  confidence  limits  for  q  obtained  on  the 
basis  of  this  assumption  remain  valid  also  when  b^,...,bk  are 

random  variables.  Clearly,  the  probability  distribution  of  tla^ 

(i  =  l,...,n)  is  the  same  as  the  distribution  of  the  number  of 
successes  in  a  sequence  of  Nb^  independent  trials,  the  prob¬ 
ability  of  success  in  a  single  trial  being  ql.  Hence 

E(Nai )  =  qlMbi  ( 88) 


(Hai)  = 


Hb.q 


(1  - 


q1) 


(89) 
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From  equations  88  and  89  we  obtain 


(?)  ■  - 


bi(l  -  q1) 


2  ft A\  _  biU  ' 

W/  Nq1 


al  a2  an 

Since  the  variates  — ,  — are  independently  distributed, 

q  q  qn 

and  since  a^  is  nearly  normally  distributed  if  N  is  not  small,  we 
can  assume  with  very  good  approximation  that  the  sum 

n  a . 

£  -T  (92) 

i-i  q 

is  normally  distributed.  We  obtain  from  equations  90  and  91 


(n  a.  \  n 

e  -t)  -y  bi- 

i=l  q1/  i=l 


1  -  a 


a2  (Z  a~i)  -  £ 

\i=l  q1/  i=l 


n  b.(l  -  q1) 


For  any  positive  a  <  1  let  Aa  be  the  value  for  which 


r  - 

'-I  ■f1’ 


dt  =  a 


The  set  of  all  values  q  for  which  the  inequality 


y"  bjO  -  q  )  "  a,  n 

&  1  1 1  ‘  *»  *  *•  / 


n  bjd  -  q‘) 


^  (95) 
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is  fulfilled  forms  a  confidence  set  for  the  unknown  probability  q 
with  confidence  coefficient  a  .  However,  formula  95  cannot  be 
used,  since  it  involves  the  unknown  quantities  b,...,b  .  Since 

ai 

— t-  converges  stochastically  to  b^  as  N  -*•  <*>,  we  change  the  stan- 

q 

a .  . 

dard  deviation  of  7! — »-  only  by  a  quantity  of  order  less  than  — 

q1  a/« 

a . 

it  v/e  replace  b^  by  -j-  .  Thus,  the  set  of  values  q  that  satisfy 

q 

the  inequalities 


/n  a.(l  *  q1)  n  In  a.(l-q‘)  (96) 

‘'WjV 


is  an  approximation  to  a  confidence  set  with  confidence 
coefficient  a  • 

Denote  by  q^  the  root  of  the  equation  in  q 

n  a  . 

=  j  ~  d0  • 

1  q3 

Then  qo  converges  stochastically  to  q  as  IU  » ,  A  considerable 

simplification  can  be  achieved  in  the  computation  of  the 
confidence  set  by  substituting  qQ  for  q  in  the  expression  of  the 

a  i 

standard  deviation  of  £ — r  .  The  error  introduced  by  this  substi- 

q1 

tution  is  small  if  11  is  large.  Making  this  substitution,  the 
inequalities  uctining  the  confidence  set  are  given  by 


•o ±  £  'A  /v 

“N/i-l  Nq“  i-1  7  '  ° 


Hence,  the  confidence  set  is  an  interval.  The  upper  end  point  of 
the  confidence  interval  is  the  root  of  the  equation 


( yt; ) 


and  the  lower  end  point  of  the  confidence  interval  is  the  root  ol 
the  equation 


(99) 


NUMERICAL  LXAMPLL 

In  all  previous  examples  it  was  assumed  that  A^  (the  number  ot 

planes  returning  with  i  hits)  was  compiled  from  such  a  large 

number  ot  observations  that  they  were  not  subject  to  sampling 

errors.  If  it  is  further  assumed  that  the  probability  q  that  a 

hit  will  down  a  plane  does  not  depend  on  the  number  oL  previous 

non-destructive  hits,  it  is  possible  to  obtain  an  exact  solution 

tor  the  probability  that  a  hit  will  down  a  plane.  Here  we 

introduce  the  possibility  that  the  A  ,...,A  are  subject  to 

o  n 

sampling  errors  but  retain  the  assumption  ot  independence.  Under 
these  less  restrictive  assumptions  we  cannot  obtain  the  exact 
solution  for  q,  but  for  any  positive  number  a  <  1  we  can  construct 
tv/o  functions  of  the  data,  called  confidence  limits,  such  that 
the  statement  that  q  lies  between  the  confidence  limits  will  be 
true  100a  percent  of  the  time  in  the  long  run.  The  confidence 
limits  are  calculated  for  a  =  .95  and  .99. 


-70- 


f 


Under  the  assumptions  of  part  I,  it  was  proved  that  no  planes 
received  more  hits  than  the  greatest  number  of  hits  observed  on  a 
returning  plane.  This  is  not  necessarily  true  when  the  possi¬ 
bility  of  sampling  error  is  introduced,  but  it  is  retained  as  an 
assumption,  since  the  error  involved  is  small. 

If  the  a^  are  subject  to  sampling  error,  and  q  is  the  true  para¬ 
meter. 


£ 


i  q 


i 

-r* 

1 


(A) 


will  be  approximately  normally  distributed  with  mean  value  1  -  aQ 

In  outlining  the  steps  necessary  to  calculate  the  confidence 
limits,  the  following  hypothetical  set  of  data  will  be  used. 

Given 


A. 


N  * 

= 

500 

ai 

a 

1 

TT 

Ao 

* 

400 

ao 

a 

O 

CO 

• 

A1 

a 

40 

al 

a 

.08 

A2 

= 

25 

a2 

= 

.05 

A3 

a 

5 

a3 

a 

.01 

A4 

= 

3 

a4 

= 

.006 

A5 

= 

2 

a5 

a 

.004 

475 

first 

step 

is 

to 

find 

A  is 


equal  to  its  mean  value,  by  finding  the  positive  root  of 


We  obtain 

. 20q5  -  . 08q4  -  .05q3  -  .Olq2  -  .006q  -  .004  =  0 

q„  =  .850  . 
o 
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The  next  step  is  to  calculate  the  standard  deviation  of 
expression  A.  This  can  be  shown  to  be  approximately  equal  to 


=  .01226  . 


n  aA 

Knowing  that  £  — r-  is  approximately  normally  distributed  with 

i=l  q1 

mean  value  1  -  a  and  the  standard  deviation  o,  we  can  determine 
o 

~ai 

the  range  in  which  2^ — -  can  be  expected  to  be  100a  percent  of 

q1 

the  time  (say  95  and  99  percent)  by  determining  A  _  and  A 
such  that 


From  the  table  or  the  areas  of  a  normal  curve,  it  is  found  that 

\g5  =  1.959964 
X>9g  =  2.575829  . 


We  can  now  calculate  the  confidence  limits  for  each  value  of  a  by 
finding  the  two  values  of  q  for  which  the  equality  sign  of  the 
following  expression  holds: 


n  a . 

T  _i  -  (1  -  a  ) 
■*— ,  1  o 

1=1  q 


*  Xaa 


It  follows  that  for  each  a,  the  confidence  limits  are  the 
positive  roots  of  the  equation 


n  a . 

Z  -r- 

1--1  q 

L  -  a  +  \  a 
o  -  Aa 

a 

Xa 

.  0122678X 

a 

1  -  a  -  X  cj 
o  na 

1 

-  a  +  X  o 
o  a 

.95 

1.559964 

.  024044 

.175956 

.224044 

.99 

2.575829 

.031600 

.168400 

.231600 

For 

a  =  .95  the 

conf idence 

limits  of  q  are 
o 

the 

positive  roots 

equat  ior. 


which  reduces  to 

.175956qb  -  .08q4 
q  =  .912, 
and  equation 


which  reduces  to 

. 224044qb  -  .  Obq4 
q  =  .801. 

similarly,  for  a  =  . 99 


+  -|  »  .175956, 
<4 


-  .U5qb  -  . 0 iq ^  -  ,0U6q 


d  r 

+  ~  =  .  224U44  , 

qJ 


-  . 0 5q b  -  .Olq2  -  ,006q 


we  have 


.004  =  0 


.004  =  0 
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. 168400q5  -  . 08q4  -  .05q3  -  . Olq2  -  . 006q  -  .004  =  0 
q  =  .935 

.231600q5  -  . 08q4  -  .05q3  -  .Olq2  -  . 006q  -  .004  =  0 
q  =  .787  . 

Summarizing  the  results  we  find  that  the  95-percent  confidence 
limits  of  q  are  .801  and  .912,  and  that  the  99-percent  confidence 
limits  are  .787  and  .935. 
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PART  VII 


MISCELLANEOUS  REMARKS1 


1 .  Factors  that  may  vary  from  combat  to  combat  but  influence  the 
probability  of  surviving  a  hit.  The  factors  that  influence  the 
probability  of  surviving  a  hit  may  be  classified  into  two  groups. 
The  first  group  contains  those  factors  that  do  not  vary  from 
combat  to  combat.  This  does  not  necessarily  mean  that  the  factor 
in  question  has  a  fixed  value  of  all  combats;  the  factor  may  be  a 
random  variable  whose  probability  distribution  does  not  vary  from 
combat  to  combat.  The  second  group  comprises  those  factors  whose 
probability  distribution  cannot  be  assumed  to  be  the  same  for  all 
combats.  To  make  predictions  as  to  the  proportions  of  planes 
that  will  be  downed  in  future  combats,  it  is  necessary  to  study 
the  dependence  of  the  probability  q  of  surviving  a  hit  on  the 
factors  in  the  second  group.  In  part  V  we  have  already  taken 
into  account  such  a  factor.  In  part  V  we  have  considered  a 
subdivision  of  the  plane  into  several  equi-vulnerability  areas 
A^#...,Ajc  and  we  expressed  the  probability  of  survival  as  a  func¬ 


tion  of  the  part  of  the  plane  that  received  the  hit.  Since  the 
probability  of  hitting  a  certain  part  of  the  plane  depends  on  the 
angle  of  attack,  this  probability  may  vary  from  combat  to  combat. 
Thus,  it  is  desirable  to  study  the  dependence  of  the  probability 
of  survival  on  the  part  of  the  plane  that  received  the  hit.  In 
addition  to  the  factors  represented  by  the  different  parts  of  the 
plane,  there  may  also  be  other  factors,  such  as  the  type  of  gun 

used  by  the  enemy,  etc.,  which  belong  to  the  second  group.  There 

are  no  theoretical  difficulties  whatsoever  in  extending  the 
theory  in  part  V  to  any  number  and  type  of  factors.  To 

illustrate  this,  let  us  assume  that  the  factors  to  be  taken  into 

account  are  the  different  Farts  Alf...,Ak  of  the  plane  and  the 

different  guns  g  ,...,g  used  by  the  enemy.  Let  q(i,j)  be  the 

probability  of  surviving  a  hit  on  part  A^  knowing  that  the  bullet 


has  been  fired  by  gun  g..  .  We  may  order  the  km  pairs  (i,j)  in  a 

sequence.  We  shall  denote  q(i,j)  by  q(u)  if  the  pair  (i,j)  is 
the  u-th  element  in  the  ordered  sequence  of  pairs.  The  problem 
of  determining  the  unknown  probabilities  q(u)  (u  =  l,...,km)  can 
be  treated  in  exactly  the  same  way  as  the  problem  discussed  in 


i'i'his  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  109  and 
AMP  memo  76.7. 


part  V  assuming  that  the  plane  consists  of  km  parts.  Any  hit  on 
part  A^  by  a  bullet  from  gun  g^  can  be  considered  as  a  hit  on 

part  A^  in  the  problem  discussed  in  part  V  where  (i,j)  is  the 

u-th  element  in  the  ordered  sequence  of  pairs. 

2.  Mon-probabilistic  interpretation  of  the  results.  It  is 
interesting  to  note  that  a  purely  arithmetic  interpretation  ol 
the  results  of  parts  I  through  V  can  be  given.  Instead  of 
defining  as  the  probability  of  surviving  the  i-th  hit  knowing 

that  the  previous  i  -  1  hits  did  not  down  the  plane,  we  define 

as  follows:  Let  l“P  be  the  number  of  planes  that  received  at  least 

i  hits  and  the  i-th  hit  did  not  down  the  plane,  and  let  M.  be  the 

total  number  of  planes  that  received  at  least  i  hits.  Then 

Mi 

=  —  .  Thus,  is  defined  in  terms  of  what  actually  hap¬ 

pened  in  the  particular  combat  under  consideration.  To  distin¬ 
guish  this  definition  of  q^  from  the  p>robabi  1  ist ic  definition,  we 

M . 

shall  denote  the  ratio  ^  by  CJ^*  The  quantity  q  is  unknown, 

since  we  do  not  know  the  distribution  of  hits  on  the  planes  that 
did  not  return.  However,  it  follows  from  the  results  of  p<art  1 
that  these  quantities  must  satisfy  equation  26.  If  we  can  assume 
that  in  the_particular  combat  under  consideration  we  have 
q^  =  ...  =  qn  then  the  common  value  q  of  these  quantities  is  the 

root  of  the  equation 


Assuming  that  q^  2.  ^2  —  *  *  *  —  '  t*lc  mlnimum  value  0°  of  fP 

derived  in  parts  III  and  IV  can  be  interpreted  as  the  minimum 
value  of  Qi  =  q^  ...  . 

The  minimum  and  maximum  values  of  Qi  derived  in  part  IV  can  also 

be  interpreted  as  minimum  and  maximum  values  of  Q^=  q^...  if 

we  assume  that  the  inequalities  X^q^  £  ^j+1  —  ^  = 

are  fulfilled.  Similarly,  a  pure  arithmetic  interpretation  of 
the  results  of  part  V  can  be  given. 
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3.  The  case  when  Y(i)  is  unknown.  In  part  V  we  have  assumed 
that  the  probabilities  y(  1 Y(k)  are  known.  Since  the 
exposed  areas  of  the  different  parts  depend  on  the 

angle  of  attack,  and  since  this  angle  may  vary  during  the  combat, 
it  may  sometimes  be  difficult  to  estimate  the  probabilities 
Y( 1 Y(k) .  Thus,  it  may  be  of  interest  to  investigate  the 
question  whether  any  inference  as  to  the  probabilities 
q(  1)  , . . .  ,q(k)  can  be  drawn  when  Y(  1 ),...,  "Y(k)  are  entirely  unknown. 
We  shall  see  that  frequently  a  useful  lower  bound  for  q(i)  can 
still  be  obtained.  In  fact,  the  value  q*(i)  of  q(i),  calculated 
under  the  assumption  that  the  parts  A ^  ( j  ji  i)  are  not  vulnerable 

(q(j)  =  1) ,  is  certainly  a  lower  bound  of  the  true  value  q(i). 
Considering  only  the  hits  on  part  A^ ,  a  lower  bound  of  q*(i),  and 

therefore  ai.so  of  q(i),  is  given  by  the  root  of  the  equation 


(100) 


where  a*  (r  =  0,1,..., n)  is  the  ratio  of  the  number  of  planes 
returned  with  exactly  r  hits  on  part  A^  to  the  total  number  of 
planes  participating  in  combat. 


The  lower  limit  obtained  from  equation  100  will  be  a  useful  one 
if  it  is  not  near  zero.  The  root  of  equation  100  will  be 

n 

considerably  above  zero  if  £  a*  is  not  very  small  as  compared 

r=i  r 

with  1  -  a*.  This  can  be  expected  to  happen  whenever  both  Y(i) 
and  q(i)  are  considerably  above  zero. 
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PART  VIII 


VULNERABILITY  OF  A  PLANE  TO  DIFFERENT  TYPES  OF  GUNS1 

In  part  V  we  discussed  the  case  Where  the  plane  is  subdivided 
into  several  equi- vulnerability  areas  (parts)  and  we  dealt  with 
the  problem  of  determining  the  vulnerability  of  each  of  these 
parts.  It  was  pointed  out  in  part  VII  that  the  method  described 
in  part  V  can  be  applied  to  the  more  general  problem  of  esti¬ 
mating  the  probability  q(i,j)  that  a  plane  will  survive  a  hit  on 
part  i  caused  by  a  bullet  fired  from  gun  j.  However,  this  method 
is  based  on  the  assumption  that  the  value  of  Y(i»j)  is  known 
where  Y(i,j)  is  the  conditional  probability  that  part  i  is  hit  by 
gun  j  knowing  that  a  hit  has  been  scored.  In  practice  it  may  be 
difficult  to  determine  the  value  of  Y(i,j)  since  the  proportions 
in  which  the  different  guns  are  used  by  the  eneihy  may  be  unknown. 
On  the  other  hand,  it  seems  likely  that  frequently  we  shall  be 
able  to  estimate  the  conditional  probability  Y( i I j )  that  part  i 
is  hit  knowing  that  a  hit  has  been  scored  by  gun  j.  The  purpose 
of  this  memorandum  is  to  investigate  the  question  whether  q(i,j) 
can  be  estimated  from  the  data  assuming  that  merely  the  quan¬ 
tities  Y ( i  I  j )  are  known  a  priori.  In  what  follows  we  shall 
restrict  ourselves  to  the  case  of  independence,  i.e.,  it  will  be 
assumed  that  the  probability  of  surviving  a  hit  does  not  depend 
on  the  non-destructive  hits  already  received. 

Let  6(i,j)  be  the  conditional  probability  that  part  i  is  hit  by 
gun  j  knowing  that  a  hit  has  been  scored  and  the  plane  survived 
the  hit.  Furthermore,  let  q  be  the  probability  that  the  plane 
survives  a  hit  (not  knowing  which  part  was  hit  and  which  gun 
scored  the  hit) .  Then,  similar  to  equation  82,  we  shall  have 


-  y( i' j  j  *  •  (101> 

Let  q( j )  be  the  probability  that  the  plane  will  survive  a  hit  by 
gun  j  (not  knowing  the  part  hit) .  Then  obviously 

q(j)  =  H  y  (ilj)q(i.j)  •  (102) 

1 

Let  6 ( i  |  j  )  be  the  conditional  probability  that  part  i  is  hit  by 
gun  j  knowing  that  a  hit  has  been  scored  by  gun  j  and  the  plane 
survived  the  hit.  Clearly 

f 


1This  part  of  "A  Method  of  Estimating  Plane  Vulnerability 
Based  on  Damage  of  Survivors"  was  published  as  SRG  memo  126  and 
AMP  memo  76.8. 
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6(  i I j  ) 


q(  j) 


(103) 


Y(i  I  j)g(i> j)  =  Y(i  I  j 
£Y(i|j)q(i,j) 


From  equation  103,  we  obtain 

q(i,j)  =  4fnDT  q(  j)  *  (104) 


The  quantity  6(i|j)  can  be  estimated  on  the  basis  of  the  observed 
hits  on  the  returning  planes.  The  best  sample  estimate  of  6(i|j) 
is  the  ratio  of  the  number  of  hits  scored  by  gun  j  on  part  i  of 
the  returning  planes  to  the  total  number  of  hits  scored  by  gun  j 
on  the  returning  planes.  Thus,  on  the  basis  of  equation  104,  the 
probability  q ( i , j )  can  be  determined  if  q ( j )  is  known. 

Now  we  shall  investigate  the  question  whether  q(j)  can  be 
estimated.  First,  we  shall  consider  the  case  when  it  is  known  a 
priori  that  a  certain  part  of  the  plane,  say  part  1,  is  not 
vulnerable.  Then  q ( i , j )  =  1  and  we  obtain  from  equation  104 


« 


1 


<5(11  j) 
Y(  1 1 }  ) 


q(  j) 


Hence , 


q(  j) 


Y(  1 1  j  ) 
6(11  j)  * 


(105) 


(106) 


Thus,  in  this  case  our  problem  is  solved.  If  no  part  of  the 
plane  can  be  assumed  to  be  invulnerable,  then  we  can  still  obtain 
upper  limits  for  q(j).  In  fact,  since  q(i,j)  <  1,  we  obtain  from 
equation  104 


q(  j) 


.  Y(i|j) 

-  TrrtiT 


(107) 


Denote  by  p(j)  the  minimum  of 
we  have 


Y(ili) 


with  respect  to 


1. 


Then 


q(j)  <  P(j)  •  (108) 

If  there  is  a  part  of  the  airplane  that  is  only  slightly 
vulnerable  (this  is  usually  the  case),  then  q(j)  will  not  be  much 
below  P(j).  Let  the  part  i^  be  the  part  of  the  plane  least 
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vulnerable  to  gun  j.  If  q(ij,j)  has  the  same  value  for  any  gun 

j,  then  q(j)  is  proportional  to  P(j).  Thus,  the  error  is  perhaps 
not  serious  if  we  assume  that  q{ j )  is  proportional  to  p(j),  i.e.. 


q ( j )  =  XP(j).  (109) 

The  proportionality  factor  X  can  be  determined  as  follows.  From 
equations  101  and  104  we  obtain 


8 (i, j ) 
Y  (  i !  j  ) 


q 


XP(  j) 


mi  j) 

Y  ( i  I  j  ) 


(110) 


Hence, 


XY  (  i  ,  j  )  =  q 


8(i, j )Y(i  I  j  ) 

17 


,])YUI] 
ill )  P  (  jf 


(111) 


Denote  X]  8(i,j)  by  S ( j )  -  Then, 
i 


From  equations  111  and  112  we  obtain 


(112) 


xYU.j)  = 


Since 


L  mij>  =  i, 

i 

we  obtain  from  equation  113 


8(  j) 


But 


A  £  £  Yd, j)  =  q  E 

j  i  j 


=  1. 


(113) 


(114) 


( 
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(115) 


Hence , 


x  = 


q£ 

j 


6(j) 

P(1T 


Since  6(j)  and  p(j)  are  known  quantities,  the  proportionality 
factor  x  can  be  obtained  from  equation  115.  The  probability  q  is 
the  root  of  the  equation 


a . 

-J.  = 


a 


o 


t 


exactly  j  hits  to  the  total  number  of  planes  participating  in 
combat . 

NUMERICAL  EXAMPLE 

In  part  V,  the  case  of  a  plane  subdivided  into  several  equi- 
vulnerability  areas  was  discussed,  and  the  vulnerability  of  each 
part  was  estimated.  The  same  method  can  be  extended  to  solve  the 
more  general  problem  of  estimating  the  probability  that  a  plane 
will  survive  a  hit  on  part  i  caused  by  a  bullet  fired  from  gun  j, 
if  assumptions  corresponding  to  those  of  part  V  are  made.  The 
first  three  of  the  four  assumptions  that  must  be  made  to  apply 
the  method  of  part  V  directly  are  identical  with  those  made  in 
part  V.  They  are : 

•  The  number  of  planes  participating  in  combat  is  large  so 
that  sampling  errors  can  be  neglected. 

•  The  probability  that  a  hit  will  not  down  the  plane  does 
not  depend  on  the  number  of  previous  non-destructive 
hits.  That  is,  q^  =  q2  =  ...  *  qQ  (say),  where  q^  is  the 

conditional  probability  that  the  i-th  hit  will  not  down 
the  plane,  knowing  that  the  plane  is  hit. 

•  The  division  of  the  plane  into  several  parts  is 
representative  of  all  planes  of  the  mission. 

The  fourth  assumption  necessary  to  apply  the  method  of  part  V 
directly  usually  cannot  be  fulfilled  in  practice.  It  is: 

•  Given  that  a  shot  has  hit  the  plane,  the  probability  that 
it  hit  a  particular  part,  and  was  fired  from  a  particular 
type  of  gun,  is  known. 
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1 


These  probabilities  depend  upon  the  proportions  in  which  differ¬ 
ent  guns  are  used  by  the  enemy.  To  overcome  this  difficulty  a 
method  that  does  not  depend  on  these  proportions  is  developed  in 
part  VIII.  The  assumptions  necessary  for  the  method  of  part  VIII 
differ  from  those  of  part  V  only  in  that  the  fourth  assumption  is 
replaced  by: 

•  Given  that  a  shot  has  hit  the  plane,  and  given  that  it 
was  fired  by  a  particular  type  of  gun,  the  probability 
that  it  hit  a  particular  part  is  known. 

The  information  necessary  to  satisfy  this  assumption  is  more 
readily  available,  and  in  the  numerical  example  that  follows  a 
simplified  method  is  suggested  for  estimating  these 
probabilities . 

The  Data 


The  numerical  example  will  be  an  analysis  of  a  set  of  hypotheti¬ 
cal  data,  which  is  based  on  an  assumed  record  of  dhmage  of  sur¬ 
viving  planes  of  a  mission  of  1,000  planes  dispatched  to  attack 
an  enemy  objective.  Of  the  1,000  planes  dispatched,  634  ( N ) 
actually  attacked  the  objective.  Thirty-two  planes  were  lost 
(L=32)  in  combat  and  the  number  of  hits  on  returning  planes  was: 

=  number  of  planes  returning  with  i  hits 


A  =  386 
o 

A2  =  120 
A  =  47 

A,  =  22 

A4  =  16 


(A) 


The  total  number  of  hits  on  all  returning  planes  is 

A,  +  2A-,  +  3A-,  +  4A.  +  5AC  = 

1  2  3  4  5 

120  +  2x47  +  3x22  +  4x16  +  5x11  =  399  . 


(B) 


/ 


These  399  hits  were  made  by  three  types  of  enemy  ammunition: 

« 
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x  Flak 

B2  20-mm  aircraft  cannon 

B. j  7.9-mm  aircraft  machine  gun 

and  the  hits  by  these  different  types  of  ammunition  were  also 
recorded  by  part  of  airplane  hit: 

Forward  fuselage 
C2  Engine 

Full  system 

C.  Remainder 
4 

The  necessary  information  from  the  record  of  damage  is  given  in 
table  7. 

TABLE  7 


NUMBER  OF  HITS 

OF  VARIOUS 

TYPES  BY 

PARTS 

Forward 

Fuel 

Total 

all 

parts 

fuselage , 

Engine, 

system. 

Remainder, 

C1 

C2 

C3 

C4 

Flak, 

17 

25 

50 

202 

294 

20-mm 

8 

7 

17 

18 

50 

cannon,  B2 

7 . 9-mm 

7 

13 

17 

18 

55 

machine 
gun,  B3 

Total  all  32  45  84  238  399 

types 

A  Method  of  Estimating  the  Probability  of  Hitting  a  Particular 


Part  Given  That  a  Shot  of  a  Particular  Ammunition  Has  Hit  the 


Plane! 

The  conditional  probability  that  a  plane  will  be  hit  on  the  i-th 
area,  knowing  that  the  hit  is  of  the  j-th  type,  must  be  deter¬ 
mined  from  other  sources  of  information  than  the  record  of 

^-Necessary  for  fourth  assumption. 
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damage.  Although  a  simplified  method  is  used  in  this  example, 
more  accurate  estimates  can  be  made  if  more  technical  data  is  at 
hand.  The  first  step  is  to  make  definite  boundaries  for  the 
areas  C^,  C^.  Next,  assume  that  each  type  of  enemy  fire 

,  toy  B3  has  an  average  angle  of  fire  6^,  63'  Finally, 

assume  that  the  probability  of  hitting  a  part  of  the  plane  from  a 
given  angle  is  equal  to  the  ratio  of  the  exposed  area  of  that 
part  from  the  given  angle  to  the  total  area  exposed  from  that 
angle. 


In  this  example  it  is  assumed  that  flak  (B^)  has  the  average 

angle  of  attack  of  45  degrees  in  front  of  and  below  the  plane, 
whereas  20-mm  cannon  and  7.9-mm  machine  gun  fire  both  hit  the 
plane  head-on  on  the  average.  The  area  is  so  bounded  that  it 

includes  areas  which,  if  hit,  will  endanger  the  pilot  and 
co-pilot.  Area  includes  the  engine  area  and  area  consists 

essentially  of  the  area  covering  the  fuel  tanks.  The  results  of 
computations,  based  on  the  above  assumptions,  are  assumed  to  be 

as  follows,  where  Y(CjjBj)1  represents  the  probability  that  a  hit 

is  on  part  C.  knowing  it  is  of  type  B.  (as  estimated  by  deter- 
—  3 

mining  the  ratio  of  the  area  of  to  the  total  area  as  viewed 
from  the  angle  8^  associated  with  ammunition  B ^ ) . 

(C) 


Y(C 

'V 

=  .058 

T<C1'B2) 

=  .143 

^(Cj 

|b3) 

=  .143 

y(c2 

1  B  ) 

=  .092 

v(c2!b2) 

■  .248 

y(C2 

IBj) 

=  .248 

t(c3 

'V 

=  .174 

T<C3|B2) 

=  .  303 

7(C3 

|B3> 

=  .  303 

y{c. 

IB.) 

=  .676 

?(c,Ib,) 

=  .  306 

Y(C. 

Is,) 

=  .  306 

J-This  notation  differs  from  the  previous  notation  of  part 
VIII.  In  the  first  part  of  part  VIII,  Y(ilj)  is  used  with  the 
understanding  that  the  first  subscript  refers  to  the  part  hit  and 
the  second  subscript  refers  to  the  type  of  bullet.  In  the 
numerical  example,  the  relationship  is  made  explicit  by  letting 
stand  for  the  i-th  part  (or  component)  and  B ^  for  the  j-th 

type  of  bullet.  The  same  device  is  used  throughout  this  example. 
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Computations  for  Method  of  Part  VIII 

Let  q(c^,Bj)  be  the  probability  of  surviving  a  hit  on  part  by 
gun  B j .  By  equation  104,  we  have 

6(C.  I B  . ) 

q(ci'Bj)  “  TTc^T  ’“V  •  (D) 

where  i(c  |B.)  is  the  probability  of  being  hit  on  part  C^, 
knowing  that  the  hit  was  scored  by  a  bullet  from  gun  B^  and  that 
the  plane  survived?  y(C^Ib^)  is  the  probability  of  being  hit  on 
part  C. ,  knowing  that  the  hit  was  scored  by  a  bullet  of  type  B.; 
and  q(Bj)  is  the  probability  that  a  plane  will  survive  a  hit  of 
type  B j ,  knowing  that  the  plane  is  hit.  This  can  be  estimated  by 
taking  the  ratio  of  the  number  of  hits  of  type  B_.  on  part  to 
the  total  number  of  hits  of  type  B^  on  returning  planes. 


Applying 

this  method 

to  the  table 

we 

obtain 

(E) 

5(C1,B1) 

*  .058 

5(C1lB2) 

— *■  • 

160 

5(C1|B3) 

=  .127 

6 (C2lB1) 

-  .085 

«(C2|B2) 

*“  • 

140 

« (G2  I  b3) 

=  .2  36 

* (c3Ib1) 

=  .170 

J(c3Ib2) 

• 

340 

mc3ib3) 

=  .309 

* (c4Ib1) 

=  .687 

J<c4Ib2) 

=  . 

360 

J(c4Ib3) 

=  .327 

The  final  quantity  required  to  calculate  q(C^,Bj)  by  equation  D 
is  q(Bj).  By  equation  109,  we  have 


q(B^)  =  APtBj) 


r(C.|B.) 

where  p  ( B  ^  )  is  the  minimum  of  — Tb)  respect  to  i. 


nun 


P<Bj)  = 


|6(C1lB^) 


Y(C2|Bj)  Y(C3|Bj)  Y ( C  4 | B  j )  | 

'  6 ( C 2  I B  j )  r  6 ( C 3  1 B  j C4 I B^ ) J 


PJB^)  =  min 

=  min 
=  .984 

p(B2)  =  min 

=  min 
=  .850 

p(B3)  =  min 

=  min 
=  .936 


.058  .092 

.058  '  .085  ' 

(l  ,  >1  , 

1.  143  .248 

(.160  '  .140  ' 

(.894  ,  >1 

1.143  .248 

1.127  '  .236  ' 

|>i  ,  >i  , 


.174  .6761 

.170  '  .687/ 

>1  ,  .984) 

.303  .306} 

.340  '  .360) 

.891  ,  . 850 J 

.303  .3061 

.309  '  .327| 

.981  ,  .936} 


(G) 


The  constant  multiplier  A  is  defined  by  equation  115 
6(B.  ) 


X=  q  5Z 


P(B.)  ' 


(H) 


where  6(Bj)  is  the  conditional  probability  that  a  hit  is  of  type 


The  determination  of  q  is  identical  with  the  procedure  of  part 
VII.  From  equation  26 


=  N  -  A 

o 

we  substitute  the  values  of  equation  a; 

248q5  -  120q4  -  47q3  -  22q2  -  16q  -  11  =  0 


(I) 


i 
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say)  . 


* 


The  root  is  .930  (  =  q  ,  say). 

o 

The  values  6(Bj)  are  obtained  directly  from  table  7  by  taking  the 
ratio  of  hits  of  type  on  returning  planes  to  the  total  number 
of  hits  on  returning  planes. 


6(8^ 

6(B2) 

6(B3) 


294 

399 

50 

399 

55 

399 


.  737 
.125 
.138 


(J) 


Substituting  the  results  of  equations  G,  I,  and  J  in  equation  H, 
we  obtain: 


6(B) 

-pTi^ 


=  93o  +  All  +  jl38  j 

J  | .984  .850  +  .936  j 


=  .93C  (1.0433) 

=  .9703 

Substituting  in  equation  F 


q(B1)  =  (.9703)  (.984)  =  .955 
q(B2)  -  (.9703)  (.850)  =  .825 
q(B3)  =  (.9703)  (.936)  =  .908 


(K) 


The  probabilities  q(C^,B^)  can  now  be  determined  from  equation  D 
by  using  the  values  given  in  equations  C,  E,  and  K. 


6  {  C .  |B.) 

q(Ci'Bj)  =  y(cJ|Bj)  q(Bj} 


q(c1,B1) 

= 

(  .058) 

( .955)/. 058 

= 

.955 

q(c2,B1) 

= 

(  .085) 

(  .955)/. 092 

= 

.  882 

q(c3,B1) 

= 

( .170) 

(  .955)/. 174 

= 

.933 

q(c4,B1) 

= 

( .687) 

( .955)/. 676 

= 

.971 

q(c1,B2) 

= 

(  .160) 

(  .  825)/. 143 

= 

.  923 

q(c2,B2) 

= 

(  .140) 

(  .825)/. 248 

= 

.466 

q(C3,B2) 

= 

(  .340) 

{ .825)/. 303 

= 

.926 

q(c4,B2) 

= 

( .360) 

(  .  825)/.  306 

— 

.971 

q(C1,B3, 

= 

( .127) 

( .908)/. 143 

= 

.  806 

q(C2,B3) 

= 

(  .236) 

( .908)/. 240 

= 

.  064 

q(C3,B3) 

= 

( .309  ) 

(  .908)/. 303 

= 

.  926 

q(c4,B3) 

= 

( .327) 

(  .908)/. 306 

= 

.970 

Comments  on  Results 


The  vulnerability  of  a  plane  to  a  hit  of  type  B ^  on  part  is 

the  probability  that  a  plane  will  be  destroyed  if  it  receives  a 
hit  of  type  on  part  .  Let  P(C^,Bj)  represent  this  vulner¬ 
ability.  The  numerical  value  of  P(C^,Bj)  is  obtained  from  the 
set  L  and  the  relationship 


P(Ci,Bj)  =  1  -  q(Ci,Bj)  (M) 


The  vulnerability  of  a  plane  to  a  hit  to  type  B^  on  part  is 
given  in  table  8. 

This  analysis  of  the  hypothetical  data  would  lead  to  the 
conclusion  that  the  plane  is  most  vulnerable  to  a  hit  on  the 
engine  area  if  the  type  of  bullet  is  not  specified,  and  is  most 
vulnerable  to  a  hit  by  a  20-mm  cannon  shell  if  the  part  hit  is 
not  specified.  The  greatest  probability  of  being  destroyed  is 
.534,  and  occurs  when  a  plane  is  hit  by  a  20-mm  cannon  shell 
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on  the  engine  t 

•  7 . 9-mm  machim 
conclusions  th< 
derived  by  the 

•  guides  for  loci 
prediction  of  1 


vulnerab: 


? 


Flak.  Bx 

«  20 -mm 

cannon,  B 

7 . 9-mm 
machine 
gun,  b3 

Vu lnerabi 1 ity 
to  hit  on 
specified  area 
when  type  of 
hit  is  un¬ 
specified* 


aThese  vulneral 
part  V,  and  asi 

is  hit,  kno' 

rtc^)  =  . 

^This  is  the  p 
hit,  When  neit 
specified. 
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on  the  engine  area.  The  next  most  vulnerable  event  is  a  hit  by  a 
7.9-mm  machine  gun  bullet  on  the  cockpit.  These,  and  other 
conclusions  that  can  be  made  from  the  table  of  vulnerabilities 
derived  by  the  method  of  analysis  of  part  VIII,  can  be  used  as 
guides  for  locating  protective  armor  and  can  be  used  to  make  a 
prediction  of  the  estimated  loss  of  a  future  mission. 


TABUS  8 

VULNERABILITY  OF  A  PLANE  TO  A  HIT  OF  A  SPECIFIED  TYPE 

ON  A  SPECIFIED  PART 


Forward 

fuselage 

Engine 

Fuel 

system 

Remainder 

Vulner¬ 
ability  to 
specified 
type  of 
hit  when 
area  is  un¬ 
specified 

Flak,  Bx 

.045 

.118 

.067 

.029 

.045 

20 -mm 

cannon,  B2 

.077 

.534 

.074 

.029 

.175 

7. 9-mm 
machine 
gun,  B3 

.194 

.136 

.074 

.030 

.092 

Vulnerability 
to  hit  on 
specified  area 
when  type  of 
hit  is  un¬ 

specified3 

.114 

.179 

.074 

.038 

.070b 

aThese  vulnerabilities  are  calculated  using  the  method  of 
part  V,  and  assuming  that  the  y(C^),  the  probability  that  part 

is  hit,  knowing  that  the  plane  is  hit,  are  as  follows: 

Y^)  «  .084  Y(C2)  -  .128  y(C3)  «  .212  Y(C4)  -  .576  . 

Whis  is  the  probability  that  a  plane  will  be  destroyed  by  a 
hit,  when  neither  the  part  hit  nor  the  type  of  bullet  is 
specified. 
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